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ABSTRACT 


A  phenomenological  and  parametric  study  was  conducted  to 
assess  the  effect  of  thermal  residual  stresses  on  the  stress- 
strain  response  of  a  discontinuous  fiber-reinforced  metal- 
matrix  composite  in  tensile  and  compressive  loading.  The 
material  chosen  for  this  investigation  was  the  SiC-whisker 
reinforced  A1  6061.  The  difference  between  composite  flow 
behavior  in  tension  and  compression,  as  well  as  the  effects  of 
volume  fraction,  fiber  aspect  ratio  and  fiber  spacing  were 
analyzed  within  the  framework  of  axisymmetric  finite-element 
models  to  determine  the  overall  constituative  response  of  the 
composite  and  to  solve  for  local  field  quantities  in  the  fiber 
and  the  matrix.  The  composite  was  modeled  as  a  periodic  array 
of  cylindrical  fibers  with  perfect  interfacial  bonding  and 
complete  fiber  alignment  with  the  tensile/compressive  axis.  It 
was  found  that  the  presence  of  residual  stresses  affected  the 
stress-strain  behavior  of  the  composite  by  influencing  the 
load  transfer  characteristics  between  the  matrix  and  the  fiber 
as  well  as  the  initiation  and  growth  of  the  plastic 
deformation  in  the  matrix.  The  results  of  the  study  indicated 
that  a  significant  strength  differential  exists  between 
tensile  and  compressive  loading  and  that  this  effect 
diminishes  with  increasing  volume  fraction.  The  results  also 
indicated  that  composite  stiffness,  yield  strength  and  work¬ 
hardening  rate  increase  with  increasing  volume  fraction  and 
fiber  aspect  ratio  and  that  the  variation  in  fiber  spacing  has 
a  relatively  small  effect  in  tension  and  a  large  effect  in 
compression  due  to  the  presence  of  residual  stresses .  Accession  For 
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I.  INTRODUCTION  AND  LITERATURE  SEARCH 


In  the  last  decade,  advancements  in  fabrication  technology 
have  brought  discontinuous  fiber  reinforced  metal  matrix 
composites  to  the  lead  in  the  search  for  the  ever  stiffer  and 
stronger  yet  lighter  and  cheaper  material.  The  addition  of 
silicon  carbide  whiskers  to  an  aluminum  6061  matrix  through  a 
powder  metallurgy  process  is  one  of  these  commercially 
practical  materials.  It  is  both  machinable  and  workable  with 
conventional  metal-forming  techniques  and  can  be  produced  in 
economically  feasible  quantities.  It  is  a  proven  composite 
with  routine  applications  that  include  automotive  engine 
components,  such  as  aluminum  pistons  with  selective 
reinforcement  and  pure  composite  connecting  rods.  Improvements 
in  mechanical  properties  over  the  unreinforced  alloy  include 
increased  elastic  modulus,  tensile  strength,  wear  resistance, 
fatigue  resistance,  and  strength-to-weight  ratio.  In  order  for 
industry  to  utilize  this  material  and  similar  metal  matrix 
composites  to  their  full  potential  detailed  knowledge  of  the 
factors  and  mechanisms  affecting  their  thermal,  electrical  and 
mechanical  properties  must  be  known. 

A.  THERMAL  RESIDUAL  STRESSES  IN  METAL-MATRIX  CC»4POSITES 

During  fabrication  of  the  metal  matrix  composite,  a 
coefficient  of  thermal  expansion  (CTE)  mismatch  between  the 


1 


matrix  and  the  reinforcing  fiber  produces  thermal  residual 
stresses.  The  exact  nature  of  these  stresses  and  their  effect 
on  the  mechanical  properties  of  the  composite  has  been  the 
subject  of  many  investigations.  Early  studies  [Refs.  1,  2, 
3,  4,  5,  6],  used  continuum  mechanics  to  investigate 
thermal  residual  stress  states  in  spherical  particle 
reinforced  MMCs,  and  in  continuous  fiber  MMCs . 

Arsenault  et  al.  [Ref.  1,  8,  9,  10]  used  electron 
microscopy  and  X-ray  diffraction  to  conclude  that  the  residual 
stresses  were  caused  by  CTE  mismatch  and  that  matrix 
strengthening  was  due  to  enhanced  dislocation  density  in  the 
matrix  material.  Additionally,  analytical  results  using  the 
shear-lag  theory  and  the  Eshelby  equivalent  inclusion  method 
were  compared  with  experimental  results  and  showed  theoretical 
methods  under-predicted  the  magnitude  of  strengthening.  To 
account  for  this  under-prediction,  Nardone  and  Prewo 
[Ref.  11]  modified  the  shear-lag  theory  to  account  for 
matrix-to-f iber  load  transfer  through  the  fiber  tip,  An 
analytical  prismatic  punching  model  was  later  developed  by 
Arsenault  and  Shi  [Ref.  12]  to  predict  the  increase  in 
dislocation  density  and  experimentally,  Derby  and  Walker 
[Ref.  13]  used  hardness  tests  to  verify  strengthening 
from  the  dislocations  generated  by  the  thermal  mismatch. 

Recently,  Arsenault  et  al.  [Ref.  14]  postulated 
that  the  observed  strength  in  SiC  whisker  reinforced  aluminum 
matrix  composites  can  be  accounted  for  by  increased 
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dislocation  density  and  reduced  subgrain  size  in  the  matrix 
due  to  the  difference  in  CTE  between  the  fiber  and  the  matrix. 
They  found  using  electron  microscopy  that  the  dislocation 
density  in  the  matrix  increases  with  fiber  volume  fraction  and 
decreases  with  increasing  particle  size.  Withers  et  al. 
[Ref.  15]  used  a  modified  Eshelby  method  in  an  effort 
to  provide  a  simple  and  reliable  field  model  to  predict 
composite  properties  such  as  elastic  modulus,  CTE,  residual 
thermal  stress,  interfacial  bond  strength  and  load  transfer  in 
short  fiber  reinforced  metal  matrix  composites. 

Computational  approaches  have  recently  gained  wide 
acclaim  for  their  flexibility  to  include  effects  such  as 
stress  concentration  zones  at  the  sharp  corners  of  the 
reinforcement,  and  determine  the  magnitudes  of  local  field 
quantities,  which  is  essential  to  the  study  of  thermal 
residual  stresses.  Finite  element  analysis  was  used  by  Dutta 
et  al.  [Refs.  16,  17]  to  model  the  deformation  zone  in 
SiC-A16061  MMC  arising  from  differential  thermal  contraction. 
In  the  latter  study  it  was  found  that  thermal  residual 
stresses  can  explain  the  differences  between  shear-lag  models 
and  experimentally  observed  data  and  can  account  for  the 
difference  between  proportional  limit  and  yield  strength  in 
discontinuous  fiber  composites. 

Povirk  et  al.  [Ref.  18]  conducted  a  parametric 
study  of  the  effects  of  fiber  spacing,  volume  fraction  and 
fiber  aspect  ratio.  The  composite  was  modeled  as  a  periodic 
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array  of  cylindrical  unit  cells  composed  of  a  perfectry 
aligned  thermo-elastic  fiber  and  thermo-plastic  matrix.  A 
temperature  history  was  imposed  on  the  cell  to  simulate  the 
post-fabrication  quench.  Results  indicated  that  transverse 
fiber  spacing  was  the  most  important  parameter  affecting  the 
distribution  and  magnitude  of  the  plastic  deformation  in  the 
matrix,  that  fiber  volume  fraction  affects  the  overall  level 
of  plastic  deformation  in  the  matrix  and  that  fiber  aspect 
ratio  has  a  negligible  effect. 

B.  STRESS-STRAIN  BEHAVIOR  OF  METAL  MATRIX  COMPOSITES 

Several  comprehensive  studies  have  been  conducted  to 
examine  the  dependence  of  composite  stress-strain  behavior  on 
material  properties  such  as  fiber  volume  fraction,  fiber 
aspect  ratio  and  fiber  spacing.  Studies  are  typically 
accomplished  using  two-dimensional  finite  element  analysis, 
but  more  accurate  account  of  the  presence  of  adjacent  fibers 
is  possible  with  three-dimensional  analysis. 

Christman  et  al.  [Ref.  19]  conducted  a  two- 
dimensional  finite  element  analysis  and  experimental  study  of 
deformation  in  SiC  whisker  reinforced  aluminum  composites 
under  tensile  loading.  The  composite  was  modeled  with  various 
periodic  arrays  of  aligned  fibers  with  cylindrical  matrix 
shells.  The  matrix  was  modeled  as  an  elastic  viscoplastic 
solid,  subject  to  isotropic  work  hardening,  and  the  fiber  was 
modeled  as  both  elastic  and  perfectly  rigid.  The  numerical 
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results  indicated  that  large  hydrostatic  stress  states  were 
developed  in  the  matrix  during  deformation  and  that  plastic 
deformation  occurred  first  at  the  fiber  tip  and  then  spread 
throughout  the  matrix.  The  results  of  the  parametric  study 
were  in  qualitative  agreement  with  experimental  results,  but 
the  model  with  the  more  realistic  constraints  significantly 
overpredicted  elastic  modulus,  yield  strength  and  tensile 
strength.  The  study  found  that  fiber  volume  fraction  was  the 
single  most  important  parameter  in  determining  composite 
strength,  and  that  the  effects  of  fiber  aspect  ratio  and 
spacing  are  strongly  dependent  on  the  extent  of  the  constraint 
imposed  on  the  matrix. 

Tvergaard  [Refs.  20,  21]  used  the  finite  element 
scheme  of  Christman  [Ref.  19],  a  staggered  (overlapping)  fiber 
arrangement,  to  account  for  matrix  shearing  between  adjacent 
fibers,  and  a  cohesion  zone  model,  to  account  for  fiber 
debonding.  Results  indicated  that  the  effect  of  matrix 
shearing  lowered  the  magnitude  of  strength  predictions  to  much 
closer  agreemient  with  experimental  data  and  that  fiber 
debonding  occurs  first  at  the  fiber  corner  followed  by  the 
fiber  top  and  sidewall  surfaces  respectively. 

Using  three-dimensional  finite  element  models.  Levy  and 
Papazian  [Ref.  22]  used  aligned  and  staggered  models 
in  longitudinal  and  transverse  tensile  loading.  The  results 
accurately  predicted  the  elastic  modulus  and  the  increase  in 
elastic  modulus  and  work-hardening  rate  with  increasing  fiber 
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volume  fraction  and  fiber  aspect  ratio.  Brockenbrough  and 
Surest  [Ref.  23]  used  three-dimensional  finite  element 
analysis  to  model  various  fiber  shapes  and  packing 
arrangements  to  study  the  effect  of  fiber  geometry  and 
distribution  on  the  plastic  deformation  characteristics  of  the 
composite.  Results  indicated  that  fiber  distribution  had  a 
pronounced  effect  on  transverse  tensile  properties  and  that 
the  evolution  of  matrix  hydrostatic  stresses  and  plasticity 
overshadowed  all  fiber  characteristics  in  longitudinal 
loading . 

To  study  localized  deformation  in  SiC-Al  composites, 
Arsenault  et  al.  [Ref.  24]  used  two  and  three- 
dimensional  finite  element  analysis  which  accounted  for  the 
presence  of  residual  stresses  to  study  clustering  effects  in 
the  fibers  at  large  fiber  volume  fractions  and  monitor  the 
behavior  of  plastic  zones  around  those  whiskers.  It  was  found 
that  plastic  deformation  is  initiated  within  the  region  of 
high  fiber  concentration  and  that  this  clustering  results  in 
a  reduction  of  composite  flow  stress,  therefore,  the  more 
uniform  the  fiber  distribution  the  stronger  the  composite. 

C.  STRENGTH  DIFFERENTIAL  EFFECT 

Residual  stresses  are  thought  to  be  the  source  of  the 
difference  in  yield  strength  and  flow  stress  between  tension 
and  compression  exhibited  in  metal  matrix  composites. 
Arsenault  and  Taya  [Ref.  9]  analytically  predicted  that 
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differences  in  yield  strength  between  tension  and  compression 
increase  with  increasing  volume  fraction.  Using  the  modified 
Eshelby  inclusion  method,  their  results  compared  well  to 
experimental  testing  of  the  short  fiber  MMC.  The  authors 
postulated  that  this  difference  was  due  to  a  tensile  residual 
stress  existing  in  the  matrix. 

Dutta  [Ref.  16]  using  finite  element  analysis,  determined 
that  the  positive  residual  shear  stresses  at  the  fiber  matrix 
interface,  which  results  from  the  thermal  mismatch,  reduce  the 
ability  of  the  matrix  to  transfer  load  to  the  fiber  in  tension 
and  enhance  it  in  compression.  Therefore,  the  composite 
compressive  strength  should  be  greater  than  the  composite 
tensile  strength .  Arsenault  and  Wu  [Ref .  25 ] 
experimentally  determined  the  magnitude  of  this  strength 
differential  effect  (SDE)  and  the  Bauschinger  Effect  (BE)  of 
various  SiC-Al  composites  and  compared  them  to  the  model  by 
Arsenault  and  Taya  [Ref.  9].  Results  indicated  that  the 
compressive  yield  strength  of  whisker  and  platelet  reinforced 
composites  is  larger  than  the  tensile  yield  strength  and  that 
the  magnitude  of  the  SDE  and  BE  increase  with  increasing 
reinforcement  volume  fraction. 

Povirk  et  al.  [Ref.  26]  conducted  a  finite  element 
analysis  of  the  effect  of  residual  stresses  on  deformation  and 
damage  mechanisms  in  SiC-Al  composites.  The  predicted  response 
with  residual  stresses  in  tension  and  compression  was  compared 
to  the  response  without  residual  stresses.  The  fiber  matrix 
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interface  was  also  modeled  in  tension  with  a  specified 
decohesion  parameter  to  characterize  the  separation  of  the 
interface.  Composite  cells  with  a  fiber  aspect  ratio  of  four 
and  a  cell  aspect  ratio  of  two  and  four  were  compared  and 
results  indicated  that  closer  end  to  end  spacing  promotes 
greater  flow  strength  in  compression  than  in  tension,  and  that 
thermal  residual  stresses  enhance  this  effect. 

D.  RESEARCH  OBJECTIVES  AND  OVERVIEW 

Past  research  efforts  have  failed  to  comprehensively 
examine  the  effect  of  thermal  residual  stresses  on  composite 
processing  variables  and  the  overall  stress-strain  response  of 
the  whisker  reinforced  metal-matrix  composite.  This  thesis, 
and  the  thesis  [Ref.  27]  that  preceded  it,  have 
attempted  to  take  an  in-depth  look,  through  finite  element 
analysis,  at  residual  stress  effects  on  the  SiC„-Al  composite. 
The  specific  objectives  of  this  thesis  include; 

1.  Study  the  phenomenology  of  the  matrix  and  fiber 
deformation  in  the  composite  with  and  without 
residual  stresses. 

2.  Study  the  effect  of  fiber  volume  fraction,  fiber 
aspect  ratio,  and  fiber  spacing  in  presence  of 
residual  stresses. 

Chapter  II  summarizes  the  assumptions  and  methods  used  in 
the  numerical  modeling  as  well  as  the  material  testing 
procedures.  Chapter  III  contains  the  results  of  this 
investigation,  and  includes  results  from  the  previous 
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investigation  [Ref.  27]  for  comparison.  Chapter  IV  summarizes 
the  relevant  conclusions  of  the  research. 
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II.  EXPERIMENTAL  APPROACH 


A.  NUMERICAL  MODELING 

The  detailed  account  of  the  construction  and  validation  of 
the  numerical  modeling  employed  in  this  thesis  is  presented  in 
preceding  thesis  work  by  Sims  [Ref.  27].  The  following  is  a 
summary  of  the  methodology  and  assumptions  used  in  this 
follow-on  study. 

1 .  Finite  Element  Modeling 

The  multipurpose  finite-element  code,  ADINA,  was  used 
to  predict  macroscopic  stress-strain  behavior,  and  localized 
microscopic  stress-strain  behavior  of  the  metal-matrix 
composite.  Composite  and  fiber/matrix  stress-strain  curves  and 
iso-stress/ iso-strain  contour  plots  were  produced.  Two- 
dimensional  axisymmetric  analysis  was  used  in  conjunction  with 
an  Updated-Lagrangian- Jaumann  (ULJ)  kinematic  formulation 
which  uses  virtual  work  to  calculate  incremental  displacements 
and  ultimately  elemental  stresses.  The  non-linear  response  of 
the  material  is  approximated  through  an  incremental  approach 
based  on  the  Newton-Raphson  equations  [Ref.  28],  that 
predicts  the  configuration  of  the  body  at  the  next  solution 
step.  The  numerical  model  used  two-dimensional  iso-parametric 
elements  and  second-order  Gaussian  integration  was  used  for 
solution  accuracy. 
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The  aluminum  matrix  was  characterized  by  a  bilinear 
stress-strain  curve,  and  was  assumed  to  undergo  isotropic  work 
hardening  and  obey  the  Von  Mises  criterion.  All  material 
properties  were  assumed  to  be  constant  within  the  temperature 
range  considered  in  the  study.  For  aluminum  6061  in  the  T6 
condition,  the  relevant  thermo-mechanical  properties  are; 
elastic  modulus  E.g,_  =  68.3  GPa  and  E„_„._  =  69.7  GPa,  work¬ 
hardening  modulus  E.=656  MPa,  yield  strength  <5^-216  MPa, 
poisson  ratio  v=0.33,  and  CTE=23 . O^K'- ;  the  SiC  fiber  values 
are  E  =  450  GPa,  v=0.2,  CTE  =  4.3  |iK'-. 

Unit  cell  configuration  is  similar  to  the  model  used 
by  Christman  et  ai.  [Ref.  19]  for  an  alig.ned  fiber  MMC  (Figure 
la)  .  The  material  is  assumed  to  be  composed  of  a  periodic 
array  of  cylindrical  unit  cells  perfectly  aligned  with  the 
tensile/compressive  axis.  Each  cell  (Figure  lb)  is  composed  of 
a  cylindrical  fiber  of  radius  R.  and  length  L,  surrounded  by 
a  matrix  shell  of  radius  R.  and  length  L^..  The  three- 
dimensional  problem  is  reduced  to  two-dimensional  by  using  an 
axisymmetric  modeling  approximation  which  assumes  that  a 
structure  with  rotational  symmetry  about  an  axis  in  the 
presence  of  a  rotationally  symmetric  load  can  be  fully 
described  using  a  two-dimensional  model  with  a  unit  radian 
width  [Ref.  29].  It  can  be  seen  from  figure  lb,  that 
the  axisymmetric  model  is  an  appropriate  analog  to  the  three- 
dimensional  array  of  hexagonal  unit  cells  [Ref.  19] . 
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Material  Representation:  Diagram  Showing  (a)  the 
Aligned  Fiber  Arrangement  and  (b)  a  Depiction  of 
the  Fiber  Imbedded  in  the  Matrix  Shell. 


2 .  Boundary  Conditions 


The  ADINA  code  dictates  displacement-based  boundary 
conditions  be  specified  along  with  the  master  degrees  of 
freedomi  for  the  entire  finite-element  model.  A  cylindrical 
coordinate  system  was  used  with  displacements  and  u_. 
corresponding  to  the  z  and  r  axes,  respectively.  For  the 
aligned  fiber  model,  the  boundary  conditions  required  for 
sym.m^etry  are  u,  =  C  along  z  =  0  and  u.  =  0  along  r  =  0.  The  boundary 
conditions  required  for  constraint  by  adjacent  unit  cells  are 
that  u.  on  z  =  L.  and  u.  on  r=R^  be  constrained  to  preserve  the 
right  circular  cylindrical  shape  of  the  unit  cell,  therefore 
the  top  and  lateral  surfaces  of  the  cell  remain  planer 
throughout  the  quench  and  the  loading  process. 

Perfect  bonding  is  assumed  along  the  fiber-matrix 
interface.  This  is  accomplished  in  ADINA  by  specifying  that 
the  fiber  and  matrix  elements  adjacent  to  the  interface  share 
the  same  boundary  nodes,  thereby  precluding  the  possibility  of 
separation  during  the  loading  sequence. 

3 .  Applied  loading 

The  loading  schem.e  used  in  the  numerical  modeling 
combines  thermal  and  displacement  loads  in  a  sequence  designed 
to  approximate  those  experienced  by  the  composite  during 
fabrication  and  subsequent  compressive  and  tensile  loading. 
The  post- fabi icat ion  quench  is  accomplished  by  subjecting  all 
nodes  in  the  unit  cell  to  a  505"C  temperature  drop  from  the 
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solutionizing  temperature  (530°C)  to  room  temperature  (25°C) 
in  a  single  solution  step,  which  eliminates  any  quench  rate  or 
bulk  thermal  effects.  The  resulting  nodal  displacements  from 
the  quench  are  stored  in  a  restart  file  to  be  used  as  the 
initial  nodal  displacements  during  the  tensile  and  compressive 
analysis . 

Tensile  and  compressive  loading  of  the  model  is 
accomplished  by  applying  an  axial  displacement  of  aL,=  (e*^!! )  L^. 
where  £  is  the  applied  true  strain,  is  post-quench  cell 
length,  and  aL^  is  applied  to  the  top  surface  of  the  unit 
cell.  The  composite  stress-strain  curves  are  constructed  by 
applying  a  aLc  corresponding  to  a  selected  value  of  £  to  the 
model.  The  average  axial  stress  was  calculated  using  a  volume 
average  approach  as  follows:  0^^,.=  (1/V.)  Jo„.dV;+ (l/V^) 
where  subscripts  c, f,  and  m  correspond  to  composite,  fiber, 
and  matrix  respectively.  In  this  equation  each  represents 
an  average  of  the  four  appropriate  Gaussian  integration  points 
within  each  element  of  the  cell.  The  average  elemental 
stresses  are  then  multiplied  by  their  corresponding  elemental 
volumes  and  summed  over  the  entire  volume  of  the  composite. 

4 .  Grid  Independence 

A  Full-Newton  iterative  scheme  was  selected  for  use  in 
this  study,  the  out-of-balance  energy,  with  a  tolerance  set  at 
10'*,  was  determined  to  be  the  basis  for  convergence.  The  mesh 
for  each  model  was  uniformly  refined  until  two  successive 
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values  of  a„.  at  £-=0  (residual  axial  stress  in  the  as-quenched 
condition)  gave  results  that  agreed  to  within  10%.  This 
parameter  was  found  to  be  the  most  sensitive  to  changes  in 
mesh  size  and  therefore  was  determined  to  be  the  most  suitable 
to  establish  grid  independence. 

5 .  Model  Implementation 

The  numerical  models  discussed  above  were  modified  for 
the  particular  variable  whose  effect  was  being  evaluated.  For 
the  models  used  to  study  fiber  volume  fraction  effects,  the 
fiber  and  the  cell  aspect  ratios  were  both  fixed  at  6  and 
numerical  results  were  carried  out  for  volume  fractions  of 
10%,  20%,  and  30%.  For  the  fiber  aspect  ratio  study,  the  fiber 
volume  fraction  was  held  at  30%;  the  fiber  aspect  and  cell 
aspect  ratios  were  set  equal  to  one  another  and  numerical 
results  were  carried  out  for  values  of  1,  3,  6,  and  10.  For 
the  fiber  spacing  study,  the  fiber  volume  fraction  was  again 
held  at  30%;  the  fiber  aspect  ratio  was  held  at  5,  and 
numerical  results  were  carried  out  for  a  cell  aspect  ratio  of 
3,  6,  and  10.  The  models  used  to  approximate  experimental  data 
were  assigned  values  of  3  for  both  fiber  cell  and  fiber  aspect 
ratios  and  numerical  results  were  carried  out  for  volume 
fractions  of  20%  and  30%. 

Composite  stress-strain  curves  and  appropriate  fiber 
loading  curves  were  produced  in  tensile  and  compressive 
loading,  with  and  without  residual  stresses.  Iso-stress 
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contour  plots  of  Von  Mises  effective  stress,  axial  stress,  ai-  i 
hydrostatic  stress  in  addition  to  iso-strain  contour  plots  of 
cumulative  effective  plastic  strain  were  generated  in  tension 
and  compression,  with  and  without  residual  stresses.  The 
aforementioned  results  were  compared  and  contrasted  to 
determine  the  effects  of  residual  stresses  on  the  metal-matrix 
composites 

B.  MA.TERIAL  TESTING 

Two  different  types  of  SiC  whisker  reinforced  composite 
material  were  tested  in  tensile  and  compressive  loading  to 
extract  mechanical  properties  for  the  use  in  the  numerical 
models  and  to  serve  as  a  comparison  to  assess  the  accuracy  of 
the  numerical  predictions.  The  composites  were  procured  from 
Advanced  Composite  Materials  Corporation  (ACMC) ,  Greer,  SC, 
and  were  produced  using  powder  metallurgy  (P/M)  techniques. 
Test  samples  were  made  by  the  manufacturer  from  a  20  vol.% 
SiC/Al  6061  MMC  with  a  extrusion  ratio  of  11:1  and  a  30  vol.% 
SiC/A16061  MMC  with  an  extrusion  ratio  of  6.5:1,  and  both 
materials  were  heat  treated  to  the  T6  condition.  Unreinforced 
A1  6061-T6  testing  samples  were  machined  from  plate  stock. 

All  specimens  were  cut  from  extruded  stock  such  that  the 
longitudinal  axis  of  the  sample  was  aligned  with  the  extrusion 
direction.  Samples  were  mounted  in  a  screw-driven  Instron 
machine  and  tested  at  a  rate  of  0.127  mm/min.  Tensile 
specimens  with  1.0"  gage  lengths  were  tested  with  clip-on 
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extensometers  to  record  engineering  strain  values.  Compression 
samples  with  1.0"  gage  lengths  and  length  to  thickness  ratios 
of  approximately  2.5,  to  avoid  buckling,  were  tested  without 
clip-on  extensometers.  A  more  accurate  means  of  testing 
materials,  to  be  compared  in  compression  and  tension,  would 
include  using  identical  test  specimens  with  extensometers  in 
a  testing  rig  designed  for  both  tension  and  compression. 

Scanning  electron  microscopy  was  used  to  obtain 
micrographs  of  the  20  and  30  vol.%  MMC,  to  determine  average 
values  of  the  fiber  aspect  ratio  and  cell  aspect  ratio,  and  to 
serve  as  a  microstructural  comparison  for  the  numerical 
modeling . 
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III.  RESULTS 


A.  EFFECTS  OF  THERMAL  RESIDUAL  STRESSES: 

A  PHENOMENOLOGICAL  STUDY 

1.  Constitutive  Composite  Response 

The  primary  aim  of  this  study  has  been  to  assess  the 
role  of  thermal  residual  stresses  (TRS)  on  the  stress-strain 
behavior  of  the  metal-matrix  composite.  The  results  contained 
in  this  section  will  describe  the  behavior  for  a  composite 
with  10  vol.%  SiC  whiskers  with  a  fiber  aspect  ratio  of  6  and 
a  composite  cell  aspect  ratio  of  6.  Figure  2a  compares  the 
aligned-f iber  model  with  and  without  thermal  residual  stresses 
(TRS)  in  both  uni-axial  tensile  and  uni-axial  compressive 
loading.  The  numerical  model  predicts  a  Young's  modulus  (E)  in 
tension  of  76.2  GPa,  a  modulus  in  compression  of  88.8  GPa,  and 
a  modulus  of  94.8  GPa  without  TRS.  The  model  predicts  a  yield 
strength  (YS)  of  368  MPa  in  tension,  385  MPa  in  compression, 
and  400  MPa  without  TRS.  Most  importantly,  the  model  shows  a 
significant  decrease  in  flow  stress  in  tension  and  a 
corresponding  increase  in  flow  stress  in  compression  when 
compared  to  the  model  without  TRS.  Additionally,  Figure  2b 
shows  that  in  both  loading  cases  with  TRS  present,  the 
proportional  limit  is  smaller  and  the  work-hardening  rate  is 
slightly  higher  when  compared  to  loading  without  TRS. 
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Figure  2.  Residual  Stress  Effects:  Effects  of  Thermal 
Residual  Stresses  (TRS)  in  Tension  and 
Compression,  (a)  Composite  Stress-Strain  Behavior 
and  (b)  Fiber  and  Matrix  Stress-Strain  Behavior. 
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When  TRS  are  not  present,  the  constitutive  stress- 
strain  behavior  of  the  composite  is  essentially  the  same  in 
both  tension  and  compression.  Figure  2b  shows  the  matrix, 
after  post-fabrication  quench,  is  initially  at  an  average 
tensile  stress  of  88  MPa;  therefore  the  start  of 
macroscopic  matrix  plastic  deformation  begins  at  a  smaller 
strain  level  in  tensile  loading  than  in  compressive  loading, 
where  the  matrix  m.ust  undergo  relaxation  from  its  tensile 
state  prior  to  a  build-up  in  compression  and  subsequent 
macroscopic  plastic  deformation  of  the  matrix.  This  may 
explain  the  smaller  Young's  modulus  observed  in  compressive 
loading  with  TRS  present. 

It  must  be  noted  that  an  initial  transient  occurs, 
predicting  an  initial  composite  modulus  higher  than  expected 
in  tensile  loading  and  lower  than  expected  in  compressive 
loading.  This  could  be  explained  by  the  tensile  residual 
stress  imposed  on  the  matrix  by  the  post-fabrication  quench. 
In  tension,  this  residual  stress  would  increase  and  allow  for 
immediate  load  transfer  to  the  fiber,  which  will  undergo  a 
stress  reversal  from  its  initial  compressive  residual  stress 
state.  In  compression,  this  initial  tensile  residual  stress  in 
the  matrix  must  undergo  a  stress  relaxation  and  reversal  prior 
to  additional  compressive  load  transfer  on  the  fiber.  Further 
study  is  required  in  this  area  at  very  small  strains  to 
substantiate  the  preceding  explanation. 
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When  comparing  Figure  2a  and  Figure  2b,  the  composite 
0.2%  offset  yield  stress  in  tensile  loading  occurs  well  after 
the  matrix  has  reached  its  apparent  yield  point  on  the  matrix 
stress-strain  curve,  both  in  the  presence  and  absence  of  TRS . 
In  compression,  the  0.2%  offset  yield  stress  occurs  prior  to 
the  matrix  reaching  its  apparent  yield  point  on  the  matrix 
stress-strain  curve.  It  should  be  noted  from  Figure  1,  that  in 
this  region  the  composite  stress-strain  curve  in  compressive 
loading  crosses  both  the  tensile  loading  curve  and  the  curve 
without  TRS,  while  there  is  no  indication  of  this  trend  in  the 
fiber  and  matrix  curves  of  Figure  2b. 

Thermal  residual  stresses  induce  a  positive  residual 
shear  stress  at  the  fiber  matrix  interface.  This  residual 
shear  stress  is  added  to  any  applied  shear  stress  from  an 
applied  load,  negative  in  tensile  loading  and  positive  in 
compressive  loading.  The  absolute  value  of  this  quantity  is 
indicative  of  the  amount  of  matrix-to-f iber  load  transfer  in 
the  composite  in  uni-axial  loading.  In  the  case  of  two  equal 
but  opposite  loadings,  it  is  clear  that  the  composite  under 
compression  should  exhibit  the  largest  load  transfer,  and 
subsequently  the  largest  strength. 

Figure  2a  illustrates  that  the  flow  stress  of  the 
composite  is  greater  in  compression  than  in  tension  throughout 
the  plastic  region  of  the  composite  response  when  residual 
stresses  are  present.  This  suggests  that  either  matrix-to- 
fiber  load  transfer  or  matrix  work-hardening  plays  a 
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significant  role  in  the  strength  difference  of  MMCs  . 
compression  and  tension.  From  Figure  2b  it  is  apparent  that 
load  transfer  may  play  the  important  role  in  the  stress-strain 
behavior  of  the  composite.  Under  tensile  loading,  the  fiber 
carries  less  load  and  the  matrix  carries  a  greater  percentage 
of  the  load  than  the  MMC  without  TRS .  Under  compressive 
loading,  this  is  reversed,  the  fiber  carries  more  load  and  the 
matrix  carries  less  load  than  the  MMC  without  TRS, 

The  work-hardening  rate  of  the  matrix  is  higher  in 
both  tension  and  compression  with  TRS  present.  The  composite 
work-hardening  rate  is  slightly  greater  in  tension  and  nearly 
the  same  in  compression  as  the  MMC  without  TRS  present,  which 
is  probably  due  to  the  initiation  of  macroscopic  plastic 
deformation  in  the  matrix  at  a  smaller  strain  in  tension  than 
in  compression. 

2.  Plastic  Zone  Initiation  and  Growth: 

Von  Mises  Effective  Stress 

Figure  3  shows  one  quadrant  of  the  unit  cell  of  a  10 
vol.%  whisker  reinforced  composite  model  with  a  fiber  aspect 
ratio  and  composite  cell  aspect  ratio  of  six  (Aj=6,A^=6),  at 
selected  levels  of  uni-axial  tensile  strain  with  and  without 
TRS.  Iso-stress  contours  of  Von  Mises  effective  stress  are 
plotted  within  the  model  cells  to  show  regions  of  the  matrix 
which  are  plastically  deforming.  The  contour  plot  of  residual 
effective  stress  (£,=0.0)  reveals  that  a  plastically  deformed 
zone  (0^;.>275  MPa)  exists  in  the  matrix  along  the  entire 
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fiber-matrix  interface  with  the  highest  stress  occurring  in 
the  vicinity  of  the  fiber  corner.  The  contour  plot  also 
indicates  plastic  deformation  has  occurred  throughout  the 
matrix  region  adjacent  to  the  lateral  fiber  sidewall.  This 
stress  is  highest  at  the  lateral  fiber-matrix  interface  and 
decreases  towards  the  cell  boundary.  The  effective  stress 
levels  in  the  elastic  fiber  are  characterized  by  a  large 
stress  zone  in  excess  of  500  MPa  in  the  middle  of  the  fiber, 
and  a  relatively  low  effective  stress  near  the  tip  of  the 
fiber.  Physically,  the  10  vol.%  MMC  has  a  large  volume  of 
matrix  material  surrounding  the  fiber.  Upon  quenching,  the 
matrix  contracts  around  the  fiber  resulting  in  a  high  residual 
compressive  stress  in  the  fiber  and  a  tensile  residual  stress 
in  the  matrix  surrounding  the  fiber.  Additionally,  as  the 
matrix  tries  to  relieve  the  residual  stresses,  it  will  deform 
plastically  at  the  vicinity  of  the  interface. 

At  a  small  tensile  strain  (£^=+0.001),  the  plastic  zone 
along  the  lateral  surface  of  the  interface  begins  to  expand  in 
the  radial  direction.  The  matrix  region  adjacent  to  the  top  of 
the  fiber  experiences  a  stress  relaxation  as  load  is 
transferred  from  the  matrix  to  the  fiber  across  the  fiber-tip 
and  end  surface.  This  effect  is  evident  in  the  growth  of  a 
high  stress  region  at  the  tip  of  the  fiber,  indicative  of  the 
fiber  bearing  a  greater  portion  of  the  applied  load  at  the 
fiber-tip.  At  moderate  strains  (£,=+0.0045),  the  effective  Von 
Mises  stress  in  the  matrix  increases  next  to  the  fiber-tip,  as 
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well  as  the  lateral  interface,  indicating  growth  in  the 
plastic  zone  size  which,  at  further  strain,  advances  to  cover 
the  entire  matrix. 

Figure  4  shows  the  same  quadrant  of  the  10  vol.%  MMC 
detailed  above,  with  iso-stress  contours  of  Von  Mises 
effective  stress  plotted  within  the  MMC  cell  model  at  selected 
levels  of  uni-axial  compression  with  and  without  TRS .  The  MMC 
without  TRS  exhibits  nearly  identical  behavior  in  both  tension 
and  compression.  Although  the  loading  is  in  the  opposite 
direction,  the  effective  stress  contour  plots  are  nearly 
identical . 

At  a  small  compressive  strain  (£^=-0.001),  the  behavior 
of  the  10  vol.%  MMC  is  different  from  that  under  tensile 
loading,  the  effective  Von  Mises  stress  next  to  the  lateral 
fiber-matrix  interface  decreases,  probably  suggesting  an 
increase  in  the  matrix-to-f iber  load  transfer  across  the  side 
wall  of  the  fiber.  The  plastically  deformed  matrix  region 
adjacent  to  the  end  of  the  fiber  begins  to  expand  towaras  the 
top  of  the  cell.  Contrary  to  the  tensile  case,  TRS  in 
compressive  loading  enhance  load  transfer  across  the  lateral 
fiber-matrix  interface,  but  do  not  result  in  enhanced  load 
transfer  at  the  fiber-tip.  At  a  moderate  compressive  strain 
(£^=-0.0045),  the  lateral  matrix  adjacent  to  the  fiber 
sidewall  continues  its  relaxation.  The  effective  stress  in 
this  region  near  the  middle  of  the  cell  is  well  beneath  the 
yield  stress  of  the  matrix  (aej;<275  MPa)  ,  and  is  indicative  of 
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with  residual  stress 


0  £c  =  -0.001 


=  -0.0045  £e  =  -0.008 


without  residual  stress 

Plastic  Zone  Development -Von  Mises  Effective 
Stress  in  Compression:  Contours  of  Effective 
Stress  for  a  10  vol.%  MMC  (Af^cs6)  With  and 
Without  Residual  Stress. 
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significant  matrix-to-f iber  load  transfer  through  the 
sidewall.  The  matrix  area  above  the  fiber-tip  has  entirely 
plastically  deformed,  and  the  very  high  stress  level  in  the 
matrix  region  of  the  fiber  corner  has  spread  along  the 
interface.  At  a  higher  compressive  strain  (ec;=-0 . 008 )  ,  the 
effective  stress  contour  plot  reveals  that  plastic  deformation 
in  the  matrix  has  increased  next  to  the  lateral  fiber-matrix 
interface  and  continued  next  to  the  fiber-tip  to  cover  the 
entire  composite  matrix. 

3.  Ctunulative  Effective  Plastic  Strain 

Figure  5  shows  one  quadrant  of  the  10  vol.%  MMC  with 
the  same  fiber  aspect  ratio  (A.-  =  6)  and  composite  cell  aspect 
ratio  (A. =6)  as  before.  Iso-strain  contour  lines  of  cumulative 
plastic  strain  are  plotted  within  the  model  cells  at  equal 
strain  values  (£-=±0.010)  in  tension  with  TRS,  compression 
with  TRS,  and  compression/tension  without  TRS.  The  contour 
plots  exhibit  higher  levels  of  plastic  strain  in  the  matrix 
along  the  lateral  fiber-matrix  interface  in  tension  and  higher 
levels  of  plastic  strain  in  the  matrix  region  above  the  fiber- 
tip  in  compression.  Both  models  e.xhibit  larger  plastic  strain, 
and  therefore  greater  worlc-hardening  when  compared  to  the 
model  without  TRS  present.  It  is  apparent  that  TRS  impede 
matrix  work-hardening  by  enhancing  matrix-to-f iber  load 
transfer  across  the  fiber-tip  in  tension  and  across  the 
lateral  fiber-matrix  interface,  which  relaxes  the  effective 
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with  residual 
stress 
tension 


with  residual  without  residual 

stress  stress 

compression  tension/compression 


Ec  =  +0.010  Ec  =  -0.010  Ec  *  ±0.010 


Figure  5.  Ciimulative  Effective  Plastic  Strain:  Effective 
Plastic  Strain  for  a  10  vol.%  MMC  with  Af=Ac=6  at 
an  Equivalent  Value  of  Strain  (Ec>=0.010)  in 
Tension  With  Residual  Stresses,  Compression  With 
Residual  Stresses,  and  Tension/Compression 
Without  Residual  Stresses. 
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stresses  in  the  adjacent  matrix  and  therefore  slows  the 
accumulation  of  plastic  strain. 

4 .  Fiber  Loading  in  Tension  and  Compression 

Figure  6  illustrates  the  axial  stress  along  half  of 
the  fiber  in  the  10  vol.%  MMC  (Aj=6,A^=6)  at  selected 
increments  of  tension  and  compression.  The  residual  axial 
fiber  stress  in  the  as-quenched  state,  ranges  from  -450 
MPA  at  the  fiber-tip  to  -1000  MPa  at  the  middle  of  the  fiber. 
In  tensile  loading,  the  axial  stress  along  the  fiber  relaxes 
from  its  residual  compressive  state,  and  subsequently  begins 
loading  in  tension.  In  compressive  loading,  the  fiber  exhibits 
further  compressive  loading,  but  at  a  lower  initial  loading 
rate  than  in  tension.  Notably,  the  magnitude  of  the  tensile 
stress  carried  by  the  fiber  is  significantly  lower  than  the 
magnitude  of  compressive  stress  carried  by  the  fiber  at  equal 
values  of  tensile  and  compressive  strain.  Additionally,  the 
incremental  increase  in  the  magnitude  of  the  stress  supported 
by  the  fiber  is  much  larger  in  tension  than  in  compression. 

Figure  7a,  a  comparison  of  axial  stress  along  the 
fiber  half  length  with  and  without  TRS,  shows  that  the  fiber 
with  TRS  supports  smaller  axial  tensile  stress  at  equivalent 
values  of  composite  strain  (e^)  .  Conversely,  in  compressive 
loading  (Figure  7b)  ,  the  fiber  with  TRS  supports  a  much 
greater  axial  compressive  stress  load  than  in  the  absence  of 
TRS  . 
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Figure  6.  Axial  Fiber  Stress  With  Residual  Stresses 
Present:  Axial  Stress  Carried  by  the  Fiber 
(A£=Ac=6) ,  Indicated  Levels  of  Composite  Strain 
at  Right,  Fiber  Tip  at  0.0,  Fiber  Middle  at  3.0. 

In  Figure  6,  the  magnitude  of  the  fiber-tip  axial 
stress  is  graphed  for  selected  strain  regimes.  In  tensile 
loading,  the  stress  carried  by  the  fiber-tip  is  greater  over 
all  strain  regimes  with  TRS  than  without.  The  opposite  is  true 
in  compressive  loading,  where  the  fiber-tip  stress  is  less 
with  TRS.  This  comparison  suggests  that  TRS  cause  the  MMC  to 
have  a  larger  load  transfer  rate  across  the  tip  of  the  fiber 
in  tension  and  a  smaller  load  transfer  rate  across  the  fiber- 
tip  in  compression. 
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Axial  Fiber  Stress  (MPa)  Axial  Fiber  Stress  (MPa) 


+0.008  w/o  TRS 

+0.008  w/TRS 
+0.004  w/o  TRS 

+0.002  w/o  TRS 
-^+0.004  w/TRS 
'^+0.001  w/o  TRS 

+0.002  w/TRS 

+0.001  w/TRS 

quench 


0.0  0.5  1.0  1.5  2.0  2.5  3.0 

Normalized  Fiber  Half  Length  (I/d) 

(a) 


(b) 

Figure  7 .  Axial  Fiber  Stress  With  and  Without  Residual 
Stresses:  Axial  Stress  Carried  by  the  Fiber  in 
(a)  Tension  and  in  (b)  Compression. 
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0.000-0.001  0.001-0.002  0.002-0.004  0.004-0.008 

True  Composite  Strain 


Figure  8 .  Incremental  Fiber  Tip  Stresses  With  and  Without 
Residual  Stresses:  Incremental  Fiber  Tip  Stress 

for  Indicated  Stress  Regimes  in 
Tension  and  Compression. 

5.  Axial  Stresses  in  the  Composite 

Figure  9  shows  a  quarter  cell  of  the  10  vol.%  MMC  at 
selected  values  of  uni-axial  tension  with  and  without  TRS .  The 
iso-stress  contours  of  axial  residual  stresses  at  £(.=0.0 
exhibit  a  fiber  under  purely  compressive  stress  (-1000<o^jf<- 
300  MPa) .  The  matrix  region  is  compressive  directly  above  the 
fiber,  and  tensile  next  to  the  lateral  fiber-matrix  interface. 
The  models  without  TRS  in  tension  and  compression,  with 
appropriate  sign  changes,  exhibited  essentially  the  same  iso- 
contour  plots. 
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At  a  small  tensile  strain  (£.=+0.001),  the  fiber,  stiil- 
in  compression,  undergoes  a  substantial  stress  relaxation  most 
notably  at  the  middle.  Concurrently,  the  matrix  experiences  a 
large  stress  increase  {O„^=+100  MPa)  in  the  region  above  the 
fiber-tip,  and  only  a  slight  stress  increase  in  the  lateral 
matrix-fiber  interface  region.  This  large  change  in  the  stress 
carried  by  the  middle  of  the  fiber,  and  the  coinciding  small 
change  in  stress  carried  by  the  matrix  region  adjacent  to  the 
fiber-matrix  interface  may  indicate  that  the  predominant  means 
of  matrix-to-f iber  load  transfer  in  the  MMC  is  through  the 
fiber  sidewall.  At  a  larger  tensile  strain  (£;.=+0 . 0045)  ,  the 
stress  in  the  matrix  above  the  fiber-tip  has  continued  to 
increase,  but  the  stress  carried  by  the  matrix  adjacent  to  the 
lateral  fiber-matrix  interface  has  actually  decreased.  Another 
possible  reason  for  this  decrease  in  stress  besides  matrix-to- 
fiber  load  transfer  is  the  outer  cell  wall  boundary  condition 
forcing  the  lateral  matrix  adjacent  to  the  fiber  to  relax  as 
the  cell  wall  moves  in  response  to  controlling  Poisson  ratio 
demands  in  the  upper  matrix  or  the  symmetry  requirements  from 
the  adjacent  cell.  When  compared  with  the  MMC  in  tension 
without  TRS,  there  was  some  evidence  of  a  smaller  net  increase 
in  axial  stress  in  the  matrix  region  above  the  fiber-tip  with 
TRS,  suggesting  greater  load  transfer  through  the  fiber-tip 
when  residual  stresses  are  present. 

Figure  10  shows  the  same  quarter  cell  of  the  MMC  in 
uni-axial  compression  with  and  without  TRS.  The  iso-stress 
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contours  of  axial  stress  exhibit  that  at  a  small  compressive 
strain  (e_=-0.001),  the  middle  of  the  fiber  exhibits  increased 
compressive  stress  relative  to  the  quenched  state.  The  matrix 
region  above  the  fiber  end  shows  an  increase  in  the 
compressive  stress  relative  to  the  as  quenched  state,  while 
the  lateral  matrix  area  remains  essentially  at  the  same  stress 
level  as  after  the  quench.  At  a  larger  compressive  strain 
(£.=-0.0045),  the  fiber  and  the  matrix  above  the  fiber  exhibit 
continued  increase  in  the  compressive  stress,  while  the  matrix 
next  to  the  lateral  fiber-matrix  interface  exhibits  a 
continued  reduction  in  the  tensile  stress  without  the  expected 
transformation  to  a  compressive  stress  state.  The  lack  of 
significant  tensile  stress  relaxation  in  the  lateral  matrix 
area  suggests  either  effective  matrix-to-fiber  load  transfer 
across  the  fiber  sidewall  or  the  outer  cell  wall  boundary 
condition  imposing  constraints  on  the  lateral  matrix  as 
discussed  previously.  When  compared  to  the  MMC  in  compression 
without  TRS,  there  is  some  evidence  of  a  larger  net  decrease 
in  axial  stress  in  the  matrix  region  above  the  fiber-tip  when 
residual  stresses  are  present,  indicating  less  matrix-to-fiber 
load  transfer  across  the  fiber-tip  when  TRS  are  present. 

€.  Hydrostatic  Stresses  in  the  Composite 

Figure  11  shows  a  quarter  cell  of  the  10  vol.%  MMC  at 
selected  values  of  uni-axial  tension  with  and  without  TRS.  The 
iso-stress  contours  of  hydrostatic  residual  stresses  (£^=0.0) 
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with  residual  stress 


without  residual  stress 

Figure  11.  Hydrostatic  Stress  in  Tension:  Contours  of 
Hydrostatic  Stress  for  a  10  vol.%  NMC  (Af=Ac=6) 
With  and  Without  Residual  Stress. 
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exhibit  a  large  compressive  hydrostatic  stress  region  (- 
250<a.y,<0.0  MPa)  in  the  matrix  above  the  fiber-tip,  that 
continues  to  the  lateral  matrix  region  at  the  fiber-matrix 
interface,  with  the  remainder  of  the  matrix  in  tensile 
hydrostatic  stress  ( 0  .  0<o,...,<250  MPa). 

Under  tensile  strain  (ec=+0 . 001&  +  0 . 0045)  ,  the 
compressive  matrix  region  above  the  fiber-tip  contracts  and 
becomes  tensile  while  the  matrix  region  adjacent  to  the 
lateral  fiber-matrix  interface  remains  in  compressive 
hydrostatic  stress,  suggesting  that  crack  nucleation  in  the 
MMC  will  occur  first  at  the  fiber-tip.  When  compared  to  the 
MMC  witn  TPS,  the  MMC  without  TRS  exhibits  tensile  hydrostatic 
stresses  throughout  most  of  the  matrix  region.  This  suggests 
that  residual  stresses  may  actually  inhibit  the  onset  of  crack 
nucleation  at  both  the  fiber-tip  and  lateral  f iber-matri”. 
interface . 

Figure  12  shows  the  same  quarter  cell  of  the  MMC  in 
uni-axial  compression  with  and  without  TRS.  The  iso-stress 
contours  of  hydrostatic  stress  exhibit  that  in  compressive 
strain  (£,=-0 . 001&-0 . 0045) ,  the  matrix  region  in  hydrostatic 
compression  adjacent  to  the  lateral  fiber-matrix  interface, 
next  to  the  middle  of  the  fiber,  contracts  and  becomes 
tensile.  This  suggests  crack  nucleation  via  debonding  of  the 
interface  is  likely  to  occur  at  the  lateral  fiber-matrix 
interface  close  to  the  middle  of  the  fiber  in  compressive 
loading.  When  compared  with  the  model  with  TRS,  the  MMC 
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with  residual  stress 


Figure  12. 


without  residual  stress 

Hydrostatic  Stress  in  Compression:  Contours  of 
Hydrostatic  Stress  for  a  10  vol.%  MMC  (Af=Acs6) 
in  Compression  With  and  Without  Residual  Stress. 


without  TRS  exhibits  a  hydrostatic  tensile  stress  zone  alcto 
the  entire  lateral  fiber-matrix  interface,  suggesting  TRS  may 
inhibit  debonding  at  fiber-matrix  interface. 

B .  PARAMETRIC  STUDY 

1 .  Fiber  Volume  Fraction  Effects 
a .  Stress-Strain  Behavior 

To  completely  study  the  influence  of  thermal 
residual  stresses  on  the  stress-strain  behavior  of  the  MMC, 
the  effect  of  increasing  the  fiber  volume  fraction  in  the 
composite  was  studied.  Three  aligned  fiber  models  of  the  SiC 
whisker  reinforced  MMC  (Figure  13a) ,  with  fiber  aspect  ratio 
and  composite  cell  aspect  ratio  equal  to  six  (A.=A,=6) ,  were 
constructed  and  numerically  modelled  in  uni-axial  compression 
and  tension,  with  and  without  thermal  residual  stresses  (TRS) . 
MMC  models  with  fiber  volume  fractions  of  20  and  30  percent 
were  compared  to  the  10  volume  percent  MMC  previously 
discussed.  It  is  seen  from  Figure  13b,  that  there  is  an 
overall  increases  in  stiffness  and  strength  with  increasing 
fiber  volume  fraction  and  that  residual  stresses  cause  a 
significant  difference  in  composite  properties  for  the  MMC  in 
compression  and  tension.  Most  notably,  the  figure  shows  that 
for  all  MMC  models,  the  elastic  modulus  in  compression  is 
substantially  less  than  the  modulus  in  tension  when  residual 
stresses  are  present. 
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True  Stress  (MPa) 


Af=Ac=6 


(a) 


True  Composite  Strain 
(b) 

Fibor  Volume  Fraction  Effects:  (a)  Schematic  of 
Model  Configurations  for  Volume  Fractions  of  10%, 
20%  and  30%  with  Af=Ac3=6  and  (b)  Composite 
Stress-Strain  Behavior  MMC  With  Residual  Stresses 
in  Tension  and  expression. 


Figure  13. 


Figure  14  shows  a  comparison  of  the  stress-strain 
behavior  of  the  10,  20  and  30  vol.%  MMC,  with  and  without 
residual  stresses.  The  10  vol.%  MMC  (Figure  14a)  exhibits 
greater  strengthening  in  compression  and  less  strengthening  in 
tension  with  residual  stresses  present  compared  to  the  10 
vol.%  MMC  without  residual  stresses.  The  20  vol.%  MMC  (Figure 
14b)  shows  an  intermediate  behavior,  where  little 
strengthening  is  observed  in  compression  or  tension  with 
residual  stresses  present,  compared  to  the  same  MMC  without 
residual  stresses.  The  30  vol.%  MMC  (Figure  14c),  exhibits 
greater  strengthening  it  tension  and  less  strengthening  in 
compression  when  residual  stresses  are  present,  which  is 
contrary  to  the  behavior  of  the  10  vol.%  MMC. 

To  understand  this  significant  fiber  volume 
fraction  effect,  the  30  vol.%  MMC  has  been  studied  and 
compared  to  the  10  vol.%  MMC  detailed  in  the  previous  section. 
The  30  vol.%  MMC  has  a  Young's  modulus  (E)  ,  when  residual 
stresses  are  present,  of  144  GPa  in  tension,  and  89  GPa  in 
compression  and  when  residual  stresses  are  absent,  a  modulus 
of  131  GPa  in  both  tension  and  compression.  The  0.2  %  offset 
yield  stress  (YS)  is  780  MPa  in  tension,  1130  MPa  in 
compression  and  720  MPa  without  TRS .  Figure  14c  shows  the  MMC 
in  compressive  loading  with  TRS,  with  its  significantly 
smaller  modulus,  increased  yield  strength  and  reduced  flow 
stress . 
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Figure  14 .  Fiber  Volume  Fraction  Effects  With  and  Without 
Residual  Stresses :  Composite  Stress-Strain  Curves 
in  Tension  and  Compression  for  MMC  with  Af3^c=6 
and  Fiber  Volume  Fraction  of  (a)  10%,  (b)  20%  and 
(c)  30% 
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Figure  15  shows  the  stress-strain  behavior  of  the 
fiber  and  matrix  of  the  10  vol.%  MMC  (Figure  15a)  and  the  30 
vol.%  MMC  (Figure  15b)  .  The  30  vol.%  MMC  fiber  carries 
approximately  the  same  stress  both  in  the  presence  and  absence 
of  TRS,  in  tension  and  compression.  This  is  very  different 
from  the  behavior  of  the  10  vol.%  MMC  where  the  fiber  in  the 
as-quenched  MMC  carries  more  stress  in  compression  and  less 
stress  in  tension  relative  to  the  case  without  TRS. 
Accordingly,  the  behavior  of  the  matrix  is  also  different  in 
the  30  vol.%  MMC.  The  matrix  in  the  as-quenched  30  vol.%  MMC 
carries  slightly  more  stress  in  tension  relative  to  the  case 
without  TRS  (similar  to  the  10  vol.%  MMC).  In  compression, 
however,  the  matrix  with  TRS  shows  a  substantially  lower 
instantaneous  modulus  (or  loading  rate,  dO-Zde^)  than  in  the 
case  without  TRS,  until  £(.*-0.0125.  This  large  difference  in 
do./de,  between  the  cases  with  and  without  TRS  in  the  30  vol.% 
MMC  and  the  absence  thereof  in  the  10  vol.%  MMC  is  one  reason 
for  the  observation  that  the  composite  flow  stress  in 
compression  is  less  than  that  in  tension  in  the  30  vol.%  MMC 
(contrary  to  what  is  observed  in  the  10  vol.%).  To  better 
understand  the  origin  of  the  low  initial  dO^/dE^  we  will  now 
investigate  plastic  deformation  in  the  matrix  of  the  30  vol.% 
MMC. 
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(b) 

Figure  15.  Fiber  Volume  Fraction  Effects  in  Fiber  and 
Matrix:  Fiber  and  Matrix  Stress-Strain  Behavior 
With  and  Without  Residual  Stresses  in  Tension  and 
Compression  for  Fiber  Volume  Fractions  of  (a)  10% 
and  (b)  30% . 
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b.  Plastic  Zone  Initiation  and  Growth 


Figure  16  shows  a  quarter  cell  of  the  30  vol.%  MMC 
with  TRS  at  selected  values  of  tension  and  compression.  Iso¬ 
stress  contours  of  Von  Mises  effective  stress  have  been 
plotted  for  comparison  with  the  10  vol.%  MMC  effective  stress 
contour  plots  (Figure  3  &  4) .  The  30  vol.%  contour  plot  of  the 
residual  effective  stress  (£,-=0.0)  is  similar  to  the  10  vol.% 
MMC  with  a  plastically  deformed  matrix  region  surrounding  the 
fiber,  but  with  reduced  effective  stress  in  the  fiber  (ae;f<225 
MPa)  . 

The  30  vol.%  MMC  in  tensile  loading  (Figure  16a) 
exhibits  effective  stress  contour  plots  that  closely  resemble 
that  of  the  10  vol.%  MMC.  In  compression,  the  effective  stress 
contour  plots  (Figure  16b)  reveal  very  different  behavior  than 
the  10  vol.%  MMC.  The  entire  matrix  region  above  the  fiber-tip 
plastically  deforms  at  a  strain  of  £^=-0.001,  in  contrast  to 
the  tensile  case,  where  the  region  does  not  completely  deform 
until  a  strain  beyond  £_=+0.0045.  Additionally,  in  the  10 
vol.%  MMC  this  matrix  region  is  not  plastically  deformed  until 
£_=  0.0045  in  tension  and  compression.  This  onset  of  large- 
scale  plastic  deformation  in  the  matrix  of  the  30  vol.%  MMC 
from  the  very  early  stages  of  compressive  loading  is  the 
reason  for  the  low  initial  loading  rate  (da„/d£(.)  and  the 
consequently  lower  flow  stress  in  the  presence  of  residual 
stresses . 
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Ec  =  0.0  Ec  =  ±0.001  Ec  =  ±0.0045 


compression  with  residual  stress 

Figure  16.  Fiber  Volume  Fraction  Effects  -  Plastic  Zone 
Development:  Von-Mises  Effective  Stress  of  30 
vol . %  MMC  With  Residual  Stresses  in  (a)  Tension 
and  (b)  Compression. 


c.  Effective  Plastic  Strain 

Figure  17  shows  the  quarter  cell  models  at  a  strain 
level  of  £,=±0.010  in  tension,  compression  and  without 
residual  stresses  for  the  10  vol.%  MMC  (Figure  17a)  and  the  30 
vol.%  MMC  (Figure  17b)  with  iso-strain  contours  of  cumulative 
effective  plastic  strain.  The  30  vol.%  MMC  contour  plots 
exhibit  an  higher  overall  cumulative  effective  plastic  strain 
throughout  the  matrix,  compared  to  the  10  vol.%  MMC.  This 
suggests  that  the  overall  worlc-hardenability  ability  of  the 
matrix  when  loaded  in  both  tension  and  compression  increases 
with  larger  volume  fraction,  due  to  the  smaller  and  therefore 
more  constrained  matrix  of  the  MMC.  Additionally,  in 
compression,  the  30  vol.%  MMC  exhibits  a  lower  level  of  worlc- 
hardening  in  the  matrix  adjacent  to  the  lateral  fiber-matrix 
interface  when  compared  to  the  MMC  without  TRS,  which  is  not 
observed  in  the  10  vol.%  MMC,  and  might  indicate  more 
effective  matrix-to-f iber  load  transfer  across  the  sidewall. 

d.  Fiber  Loading 

Figure  18  shows  the  average  axial  fiber  stress 
along  half  of  the  fiber  length  in  the  10  vol.%  MMC  (Figure 
18a)  and  the  30  vol.%  MMC  (Figure  18b)  at  selected  increments 
of  tension  and  compression.  The  residual  axial  fiber  stresses 
(Ojj.)  are  significantly  smaller  in  the  30  vol.%  MMC  than  in 
the  10  vol.%  MMC.  Further,  unlike  the  10  vol.%  MMC  where  the 
magnitude  of  the  tensile  fiber  stress  supported  at  a  given 
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Ee  =  ±0.010 


Figure  17.  Cumulative  Effective  Plastic  Strain:  Effective 
Plastic  Strain  at  Equivalent  Values  of  Composite 
Strain  (e^^O.OlO)  in  Tension  with  TRS, 
Compression  With  TRS,  and  Tension/Compression 
Without  TRS  for  (a)  10  vol.%  MMC  and  (b)  30 

vol.%  MMC  with  Af^c^e. 
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riguro  18. 


Axial  Fiber  Stress  Hith  Residual  Stresses:  Axial 
Stress  Carried  by  the  Fiber  at  Indicated 
Cmnposite  Strains  at  Right  in  (a)  10  vol.%  MMC 
and  (b)  30  Vol.%  MMC. 


tensile  strain  is  less  than  the  magnitude  of  the  compressive 
fiber  stress  supported  at  the  same  compressive  strain,  the 
magnitude  of  the  tensile  stress  supported  by  the  fiber  of  the 
30  vol.%  MMC  at  a  given  tensile  strain  is  about  equal  to  the 
magnitude  of  the  compressive  fiber  stress  supported  at  the 
same  compressive  strain.  This  indicates  a  smaller  disparity  in 
the  load  bearing  ability  of  the  fiber  in  tension  and 
compression  in  the  MMC  with  increasing  fiber  volume  fraction. 

In  Figure  19,  a  comparison  of  axial  fiber  stress 
along  half  of  the  fiber  length  with  and  without  TRS  in  tension 
(Figure  19a)  and  in  compression  (Figure  19b)  shows  that  at  a 
given  strain  the  axial  fiber  stress  supported  by  the  fiber  is 
nearly  the  same  with  and  without  TRS.  This  is  not  observed  in 
the  10  vol.%  MMC,  which  in  tension  and  compression,  exhibited 
a  large  stress  between  the  presence  and  absence  of  residual 
stresses.  It  is  suggested  that  TRS  plays  a  less  important  role 
on  fiber  load-bearing  ability  with  increasing  fiber  volume 
fraction  in  the  MMC. 

Figure  20  shows,  specifically,  the  fiber-tip  axial 
stress  graphed  for  selected  strain  regimes  in  tension  and 
compression,  with  and  without  TRS.  The  ordinate  of  the  plots 
shows  the  increase  in  fiber-tip  axial  stress  (O^J  for  given 
increases  in  the  applied  axial  composite  strain  (e^,)  .  It  is 
observed  that  the  10  vol.%  MMC  (Figure  20a)  is  nearly 
identical  to  the  30  vol.%  MMC  suggesting  that  although  TRS 
significantly  alter  fiber-tip  load  transfer  in  compression  and 
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Figure  19.  Axial  Fiber  Stress  30  vol.%  MMC  With  and  Without 
Residual  Stresses:  Axial  Stress  Carried  by  the 
Fiber  in  (a)  Tension  and  (b)  Compression. 
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Figuro  20.  Incremental  Fiber  Tip  Stresses  With  and  Without 
Residual  Stress:  Incremental  Fiber-Tip  Stress 
(ozZ(fi^_(i-,/Aee)  for  Indicated  Stress  Regimes  in 
(a)  10  voi.%  MNC  and  (b)  30  vol.%  MMC. 
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tension  this  phenomenon  is  not  affected  by  fiber  volume 
fraction . 

e.  Hydrostatic  Stress 

Figure  21  shows  a  quarter  cell  of  the  30  vol.%  MMC 
in  tension  (Figure  21a)  and  in  compression  (Figure  21b)  TRS . 
The  contour  plots  of  hydrostatic  residual  stress  (e;.=  0.0) 
shows  that  the  matrix  region  above  the  fiber-tip  is  tensile 
with  the  highest  stress  at  the  fiber-matrix  interface.  This  is 
not  observed  in  the  10  vol.%  MMC,  which  has  a  large 
compressive  matrix  above  the  fiber-tip.  Under  tensile  strain 
(e.=+0 . 001&+0 . 0045) ,  both  the  tensile  matrix  region  above  the 
fiber-tip  and  the  compressive  matrix  region  adjacent  to  the 
lateral  fiber-matrix  interface  expand,  suggesting  that  crack 
nucleation  in  a  MMC  occurs  first  at  the  fiber-tip  in  tension 
as  in  the  10  vol.%  MMC. 

Under  compressive  strain  (£„=-0 . 001&-0 . 0045) ,  the 
iso-stress  contours  of  hydrostatic  stress  show  that  the 
tensile  matrix  region  above  the  fiber-tip  becomes  compressive 
similar  to  the  10  vol.%  MMC  with  a  compressive  matrix  above 
the  fiber-tip  and  tensile  matrix  region  adjacent  to  the 
sidewall  of  the  fiber,  suggesting  that  interfacial  decohesion 
is  likely  to  occur  at  the  lateral  fiber-matrix  interface  in 
compressive  loading.  The  crack  nucleation  behavior  in  the  MMC 
does  not  exhibit  any  significant  dependence  on  the  fiber 
volume  fraction. 
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±  0.0045 


2 .  Fiber  Aspect  Ratio  Effects 

a .  Stress-Strain  Behavior 

To  understand  the  relationship  between  composite 
strength  and  fiber  aspect  ratio  in  the  presence  of  residual 
stresses,  four  models  of  the  30  vol.%  MMC  were  constructed  and 
numerically  modelled  in  tension  and  compression  with  residual 
stresses  present  (Figure  22a) .  The  fiber  aspect  ratio  and  cell 
aspect  ratio  were  set  equal  to  one  another  to  eliminate 
spacing  effects  (A--=A.)  .  Aspect  ratio  values  of  one,  three, 
six,  and  10  were  modeled  to  simulate  actual  values.  Figure  22b 
shows  that  the  MMCs  exhibited  greater  stiffness,  yield 
strength,  flow  stress  and  work-hardening  rates  with  increasing 
fiber  aspect  ratio  in  tension  and  compression.  Additionally, 
it  is  observed  that  the  strengthening  is  less  pronounced  as 
the  fiber  aspect  ratio  increases  in  tension  and  compression 
also.  Comparing  the  compressive  and  tensile  stress-strain 
behavior  of  equivalent  MMC  reveals  strength  differential 
effects,  specifically  in  elastic  modulus,  similar  to  the 
effect  discussed  in  the  previous  section. 

b.  Fiber  Loading  Rate  Behavior 

Figure  23  shows  the  fiber  stress-strain  behavior 
for  10  vol.%  and  30  vol.%  MMC  with  two  equivalent  fiber-cell 
aspect  ratios  (A.=A(.)  of  three  and  six.  It  is  observed  that  the 
stress  carried  by  the  fiber  increases  with  increasing  fiber 
aspect  ratio  and  larger  fiber  volume  fraction.  In  tension,  the 
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Figure  22.  Fiber  Aspect  Ratio  Effects:  (a)  Schematic  of 
Model  Configuration  for  30  vol . %  MMC  with 
AfsAc=10,  6,  3,  and  1,  and  (b)  Composite  Stress- 
Strain  Behavior  in  Tension  and  Compression  in  the 
Presence  of  Residual  Stresses. 


fiber  loading  is  initially  dependent  on  fiber  volume  fracti  '.'.  • 
and  only  at  strains  e. >+0.004  does  the  fiber  aspect  ratio  play 
an  important  role  in  the  fiber  loading  behavior.  In 
compression,  the  fiber  loading  is  initially  dependent  on  both 
fiber  volume  fraction  and  fiber  aspect  ratio,  even  though  the 
compressive  stress-strain  behavior  of  the  fibers  at  larger 
strains  £.>-0.010  are  similar  to  the  tensile  behavior. 


Figure  23.  Fiber  Load-Bearing  Behavior:  Effect  of  Fiber 
Aspect  Ratio  and  Volume  Fraction  on  the  Stress 
Carried  by  the  Fiber  with  Residual  Stresses. 


3.  Cell  Aspect  Ratio  Effects 

The  effect  of  fiber  spacing  was  studied  using  three  30 
vol.%  MMC  models  (Figure  24a)  with  a  fiber  aspect  ratio  of 
five  (A.-=5)  and  a  composite  cell  aspect  ratios  (A^.)  of  three. 
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six,  and  10.  Figure  24b  shows  the  composite  stress-strain 
behavior  in  tension  and  compression  with  residual  stresses 
present.  In  tension,  fiber  spacing  has  no  appreciable  effect 
on  composite  stiffness.  However,  in  MMCs  with  close  end-to-end 
and  side-to-side  spacing  an  increase  in  composite  work¬ 
hardening  is  observed  due  to  additional  matrix  constraint, 
which  results  in  somewhat  larger  plastic  strains,  and 
consequently  slightly  higher  matrix  flow  stress.  In 
compression,  fiber  spacing  has  a  large  effect  on  both  the 
stiffness  and  flow  stress  of  the  MMC.  It  is  suggested  that  in 
compression  there  is  enhanced  load  transfer  through  the 
sidewall,  therefore  the  matrix  next  to  the  lateral  fiber- 
matrix  interface,  which  is  volumetrically  larger  for  A„=3 
bears  little  load,  with  the  fiber  bearing  a  correspondingly 
large  load.  Hence  A-=3  shows  the  largest  strength  in 
compression.  Likewise,  A.=10  which  has  the  least  matrix  volume 
next  to  the  lateral  interface,  should  have  the  lowest  strength 
in  compression. 

C .  EXPERIMENTAL  RESULTS 

In  order  to  numerically  model  a  composite  material 
accurately,  the  model  must  be  similar  both  geometrically  and 
materially  to  the  actual  composite.  Figure  25  shows  a  SEM 
micrograph  of  a  30  vol.  SiC„/Al  6061  composite.  The  picture 
exhibits  three  important  aspects  of  the  actuals  composite; 
First,  there  is  a  significant  degree  of  randomness  in  fiber 
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Figure  25.  SEM  Micrograph:  Scanning  Electron  Microscopy 
Photography  at  3/45  KX  of  a  30  vol.%  SiC  Fiber- 
Reinforced  A1  6061  Composite  With  a  6.5:1 
Reduction  Ratio. 


alignment.  Second,  there  is  a  non-uniformity  (clustering)  in 
fiber  distribution  throughout  the  material,  and  third,  there 
is  a  wide  range  of  fiber  aspect  ratios  encountered  in  the 
miaterial.  Takao  al.  [Ref.  29]  found  that  fiber  alignment 
had  a  greater  effect  on  composite  elastic  modulus  than  volume 
fraction.  Christman  et  al.  [Ref.  19]  found  that  clustering 
greatly  reduced  the  resistance  to  plastic  flow  in  the 
composite  by  decreasing  the  magnitude  of  matrix  triaxiality. 
Takao  and  Taya  [Ref.  30]  found  that  the  wider  the 
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range  of  aspect  ratio,  the  greater  the  error  in  the  use  of 
average  aspect  ratio  in  analytical  results. 

From  Figure  25  the  average  composite  fiber  aspect  ratio 
was  estimated  to  be  about  2.2  with  a  low  of  1.0  and  a  high  of 
22.9.  This  estimate  is  expected  to  be  low  due  to  fiber 
misalignment.  The  numerical  model  with  fiber  aspect  ratio  and 
cell  aspect  ratio  set  at  three  was  determined  to  best 
represent  the  actual  composite. 

Figure  26a  shows  the  analytical  results  in  tension  and  the 
tensile  experimental  results  for  both  the  20  and  30  vol.%  MMC, 
it  was  found  that  the  aligned  fiber  model  overpredicts  both 
the  yield  strength  and  work-hardening  rate  of  the  actual 
composite.  This  can  be  attributed  to  the  fact  that  the  model 
assumes  perfectly  aligned,  parallel  fibers,  which  is  clearly 
not  the  actual  case  (Figure  25)  .  Additionally,  the  model 
ignores  effects  of  clustering,  which  is  known  to  reduce  the 
resistance  to  plastic  flow  [Ref.  19],  and  the  interfacial 
decohesion,  which  is  probably  of  importance  in  actual 
composites . 

The  tensile  and  compressive  results  for  both  the  numerical 
models  and  the  experimental  testing  are  shown  in  Figure  26b, 
where  composite  flow  stress  versus  plastic  strain  is  plotted. 

The  numerical  model  predicts  the  higher  yield  strength  in  « 

compression  than  in  tension  and  a  higher  yield  strength  and 
work-hardening  rate  with  increasing  fiber  volume  fraction, 
similar  to  the  actual  composite  behavior.  The  20  vol.%  model 
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Figure  26.  Experimental  Testing  Results:  Comparison  of 
Experimental  Testing  and  Numerical  Results  in  (a) 
Tensile  Loading  and  (b)  Tensile  and  Compressive 
Loading  With  Plastic  Strain  on  the  Horizontal 
Axis . 
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exhibits  better  experimental  correlation  than  the  30  vol,% 
MMC,  which  can  be  expected  due  to  the  mismatch  between  fiber 
aspect  ratio  and  fiber  spacing  assumed  in  the  model  and  that 
occurring  in  the  actual  composite.  The  greater  level  of  error 
in  models  with  higher  volum.e  fractions  can  probably  be 
attributed  to  the  reduced  lateral  and  longitudinal  spacing 
between  the  fiber  and  the  cell,  which  results  in  increased 
levels  of  matrix  constraint,  and  hence  has  a  significant 
effect  on  the  plastic  flow  behavior  of  the  material. 
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IV.  CONCLUSIONS 


The  uniaxial  stress-strain  response  of  a  SiC„/Al  6061 
composite  was  investigated  to  assess  the  role  of  thermal 
residual  stresses  on  the  constitutive  behavior  of  the 
material.  The  following  summarizes  the  relevant  findings  of 
the  phenomenological  study  and  the  parametric  study: 

A.  EFFECT  OF  THERMAL  RESIDUAL  STRESSES 


1.  Elastic-plastic  and  thermo-elastic-plastic 
analysis  using  finite-element  models  of  fiber 
reinforced  metal-matrix  composites  can  be  used  to 
predict  trends  in  the  uniaxial  stress-strain 
behavior  given  the  thermo-mechanical  properties  of 
the  constituent  phases  as  a  starting  point. 
However,  exact  correlation  with  actual 
experimental  results  is  difficult  since  several 
idealistic  assumptions  which  are  made  for  modeling 
expediency  are  often  not  satisfied  in  real 
composites.  This  usually  translates  into  over¬ 
estimation  of  the  composite  strength  and  stiffness 
properties  resulting  from  the  model. 

2.  Residual  stresses  significantly  impact  the  stress- 
strain  behavior  of  the  composite,  principally  by 
affecting  the  load  transfer  characteristics 
between  the  matrix  and  the  fiber, 

3.  Residual  stresses  affect  the  initiation  and 
continuation  of  plastic  deformation  in  the  matrix 
by  producing  an  initial  plastic  zone  adjacent  to 
the  fiber/matrix  interface  that  is  affected 
somewhat  differently  in  compression  than  it  is  in 
tension . 

4.  Residual  stresses  affect  the  fiber  load  transfer 
characteristics  at  the  fiber  tip,  increasing  load 
transfer  through  the  fiber  ends  in  tension  and 
decreasing  it  in  compression. 
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5. 


Residual  stresses  resulting  in  a  negative 
hydrostatic  stress  state  in  certain  locales  of  the 
matrix  adjacent  to  the  fiber-matrix  interface  may 
impede  decohesion  at  the  interface. 


B .  PARAMETRIC  STUDY 


1.  The  stiffness,  yield  strength  and  work-hardening 
rate  of  the  composite  increase  with  increasing 
fiber  volume  fraction  and  fiber  aspect  ratio  in 
tension  and  compression. 

2.  The  strength  differential  effect  produced  by 
residual  stresses  is  less  pronounced  with 
increasing  fiber  volumie  fraction. 

3.  Variations  in  lateral  fiber  spacing  have  little 
effect  in  tension,  but  have  significant  effect  in 
compression  on  composite  properties. 
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APPENDIX  A;  SAMPLE  ADINA  INPUT  FILE 


*  THIS  INPUT  FILE  SIMULATES  THERMAL  AND  TENSILE  LOADING  OF  A 

*  METAL-MATRIX  COMPOSITE.  THE  PROGRAM  INCLLTDES  COMMENTS  AT 

*  APPROPRIATE  LOCATIONS  TO  ILLUSTRATE  SOME  OF  THE  IMPORTANT 

*  FEATURES  OF  ADINA  IN  NON-LINEAR  ANALYSIS. 

* 

*  THE  FOLLOWING  BLOCK  OF  CODE  CONTAINS  ADINA  CONTROL  INFORMATION 
DATABASE  CREATE 

CONTROL  ORIGIN=UPPERLEFT  PLOTUNIT=CM 
FILEUNITS  LIST=6  L0G=6  ECHO=6 

★ 

*  THE  FOLLOWING  BLOCK  OF  CODE  CONTAINS  SPECIFICATIONS  FOR  MODEL  DEGREES  OF 

*  FREEDOM,  SOLUTION  TIME  PERIOD  LENGTH  AND  TIME-STEP  SIZE.  THE  FIRST 

*  "MASTER"  COMMAND  LINE  IS  USED  IN  ELASTOPLASTIC  ANALYSIS;  THE  SECONT) 

*  COMMAND  SIMULATES  THE  QUENCH  PROCESS  AND  THE  THIRD  SIMULATES  THE 

*  TENSILE  LOADING  FOLLOWING  THE  QUENCH. 

*  CAN  BE  USED  IN  DURING  A  GIVEN  RUN,  THE  OTHER  MUST  BE  COMMENTED  OUT  USING 

*  A 

MASTER  ID0F=100111  ITP56=0  NSTEP=6  DT=1.  0 
*MASTER  IDOF=100111  ITP56=2  NSTEP=1  DT=1. 0 
*MASTER  ID0F=100111  ITP56=2  NSTEP=6  DT=1. 0 

* 

*  THE  FOLLOWING  BLOCK  OF  CODE  CONTAINS  THE  KINEMATIC  ANT)  ITERATION  METHOD 

*  SPECIFICATIONS.  AUTOMATIC-ATS  IS  A  FEATURE  USED  IN  NON-LINEAR  ANALYSIS 

*  TO  ASSIST  IN  ACHIEVING  THE  REQUIRED  TOLERANCE  FOR  THE  OUT-OF-BALANCE 

*  ENTRGY  NORM. 

KINE.MATICS  DISP=LARGE  STRAINS=LARGE 
ITERATION  METHOD=FULL-NE\»TON  LINE-SEARCH=YES 
AUTOMATIC-ATS  N=10 

iv 


*  THE  FOLLOWING  COMMANDS  RELATE  TO  THE  PROGRAM  OUTPUT.  COMME.NT  IN  THE 

*  "WORKSTATION"  COM.MAND  IF  A  HARDCOPY  OF  THE  MODEL  CONFIGURATION  IS 

*  DESIRED. 

HEADING  ’VF=30%  AF=3  AC=3' 

PRINTOUT  VOLUME=MAX 

PORTHOLE  FORMATTED=YES  FILE=60 

*WORKSTATION  SYSTE.M=4  DEVICE=0  OPTION=0 

* 


*  THE  FOLLOWING  CODE  SPECIFIES 

*  NUMBERS  TO  THE  LEFT  ARE  TIME 

*  VALUES  OF  TEMPERATURE  CHANGE 
TIMEFUNCTION  1 

0.  0  0. 0 

1.0  -505.0 

* 


THE  NATURE  OF  THE  APPLIED  TEMPERATURE  LOAD. 
VALUES  AND  THE  NUMBERS  TO  THE  RIGHT  ARE  THE 
FROM  THE  INITIAL  TEMPERATURE  OF  THE  MODEL. 


*  THE  FOLLOWING  TIME  FUNCTIONS  REPRESENT  THE  DISPLACEMENT  LOADING  SEQUF:NCE 

*  THE  NUMBERS  TO  THE  RIGHT  ARE  DISPLACEMENTS  FROM  THE  REFERENCE  LENGTH  OF 

*  THE  MODEL  THAT  CORRESPOND  TO  A  GIVEN  LOGARITHMIC  STRAIN  TO  BE  USED  IN 

*  THE  STRESS-STRAIN  CURVE.  NOTE:  THE  DISPLACEMENTS  FOR  THE  MODEL  WITH 

*  RESIDUAL  STRESSES  START  AT  TIME=2. 0  SINCE  THE  QUENCH  PROCESS  OCCURRED 

*  FROM  TIME  0  to  1.0.  IN  ADDITION,  THE  DISPLACEMENTS  ARE  INITIALLY 

*  NEGATIVE  DUE  TO  THE  QUENCH  PROCESS.  THE  SUBSEQUENT  DISPLACEMENTS 
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*  ARE  ADDED  TO  THE  NEW  REFERENCE.  THE  DISPLACEMENT  SPECIFICATIONS 

*  IN  ADINA  ARE  ALWAYS  ABSOLUTE  DISPLACE.MENTS ,  MEANING  THAT  A  RELATIVE 

*  CHANGE  IN  DISPLACEMENT  CANNOT  BE  INPUT. 

* 

***  ALIGNED  TENSION  WITH  RESIDUAL  STRESS  *** 

TIMEFUNCTION  2 


0.0 

0.0 

2.0 

-.0244957 

3.0 

-.0200380 

4.0 

-.0111091 

5.0 

-.00216242 

6.0 

0.00680222 

7.0 

0.0157848 

***  ALIGNED  COMPRESSION  WITH  RESIDUAL  STRESS  *** 
TIMEFUNCTION  3 
0.0  0.0 
2.0  -.0333978 

3.0  -.0378422 

4.0  -.0467167 

5.0  -.0555754 

6.0  -.0644154 

7.0  -.0732238 

* 

***  ALIGNED  WITHOUT  RESIDUAL  STRESS  *** 
TIMEFUNCTION  4 
0.0  0.0 
1.0  0.00448224 

2.0  0.00896897 

3.0  0.0179559 

4.0  0.0269608 

5.0  0.0359837 

6.0  0.0450247 

7.0  0.0677065 


*  THE  FOLLOWING  BLOCK  OF  CODE  INPUTS  THE  NODAL  COORDINATES  AND  SPECIFIES 

*  THE  NUMBER  OF  ELEMENTS  DESIRED  BETWEEN  NODES. 

SYSTEM  1  TYPE=CARTESIAN 

COORDINATES 

ENTRIES  NODE  Y  Z 


1 

1.  49 

4.  48 

2 

1.  0 

4.  48 

3 

0.  0 

4.  48 

4 

1.  49 

3.  0 

5 

1.  0 

3.  0 

6 

0.  0 

3.  0 

7 

1.  49 

2.  24 

8 

1.  0 

2.  24 

9 

0.  0 

2.  24 

10 

1.  49 

0.  0 

11 

1.  0 

0.  0 

12 

0.  0 

0.  0 
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★ 

LINE 

STRAIGHT 

1 

2 

EL=3 

LINE 

STRAIGHT 

2 

3 

EL=6 

LINE 

STRAIGHT 

4 

5 

EL=3 

LIKE 

STRAIGHT 

5 

6 

EL=6 

LINE 

STRAIGHT 

7 

8 

EL=3 

LINE 

STRAIGHT 

8 

9 

EL=6 

LINE 

STRAIGHT 

10  11  EL=3 

LINE 

STRAIGHT 

11  12  EL=6 

*• 

LINE 

STRAIGHT 

1 

4 

EL=9 

LINE 

STRAIGHT 

4 

7 

EL=5 

LINE 

STRAIGHT 

7 

10 

EL=14 

LINE 

STRAIGHT 

2 

5 

EL=9 

LINE 

STRAIGHT 

5 

8 

EL=5 

LINE 

STRAIGHT 

8 

11 

EL=14 

LINT 

STRAIGHT 

3 

6 

EL=9 

LINE 

STRAIGHT 

6 

9 

EL=5 

LINE 

STRAIGHT 

9 

12 

EL=14 

*  THE  FOLLOWING  CODE  DEFINES  MATERIAL  PROPERTIES  TO  BE  USED  IN  THE  MODEL. 

*  THE  FIRST  MATERIAL  SPECS  CORRESPOND  TO  THE  FIBER  AND  THE  SECOND  TO  THE 

*  ALUMINUM  MATRIX.  EACH  SET  CONTAINS  THE  CORRESPONDING  COMMANDS  FOR 

*  ELASTOPLASTIC  A.NALySIS.  THE  "TREF"  COMMAND  MUST  BE  SET  TO  530.  0  FOR  THE 

*  QUENCH  RUN  AND  25.  0  FOR  THE  SUBSEQUENT  TENSILE  DISPLACEMENT  RUN. 

MATERIAL  1  THERMO-ELASTIC  TREF=530.  0 

0.0  4S0E9  0.2  4.  3E-6 

1000.0  450E9  0.2  4.  3E-6 

♦MATERIAL  1  ELASTIC  E=450E9  NU=0.  2 

* 

MATERIAL  2  THERMO-PLASTIC  TREF=530.  0 

0.0  68.3E9  0.33  276E6  656E6  23.  OE-6 

1000.0  68.  3F.9  0.33  276E6  656E6  23.  OE-6 

♦MATERIAL  2  PLASTIC  E=68. 3E9  NU=0. 33  y=276E6  ET=656E6 

* 

*  THE  FOLLOWING  CODE  SELECTS  GROUPS  OF  ELEMENTS  IN  THE  MODEL  AND  ASSIGNS 

*  THEM  PREVIOUSLY  DEFINED  MATERIAL  PROPERTIES 
EGROUP  1  TU'ODSOLID  SUB=0  M=1 

GSURF  5698  EL1=6  EL2=5  N0DES=4 

GSURF  89  12  11  EL1=6  EL2=14  N0DES=4 

* 


EGROUP  2  TWODSOLID  SUB=0  M=2 


GSURF  1254 
GSURF  2365 
GSURF  4567 
GSURF  7  8  11  ,10 

* 


EL1=3  EL2=9 
EL1=6  EL2=9 
EL1=3  EL2=5 
EL1=3  EL2=14 


N0DES=4 

N0DES=4 

N0DES=4 

N0DES=4 


♦  THE  FOLLOWING  CODE  CONTAINS  THE  BOUNDARY  CONDITION  SPECIFICATIONS. 

*  FOLLOWING  CONDITIONS  ARE  APPLICABLE  FOR  THE  ALIGNED  FIBER  MODEL. 

*  APPENDIX  B  CONTAINS  A  SECTION  OF  CODE  ILLUSTRATING  THE 

♦  APPROPRIATE  MODIFICATIONS  FOR  THE  STAGGERED  FIBER  MODEL. 


THE 
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BO  111111  NODES 
12 

BO  101111  NODES 
10  /  11 
34  TO  40 
BO  110111  NODES 
3/6/9 
89  TO  112 

★ 

CONSTRAINT 

1333  TO  2333 
13  3  3  3  TO  19  3  3  3 
1  2  10  2 
4  2  10  2 
7  2  10  2 

41  2  10  2  TO  64  2  10  2 

* 

*  THE  FOLLOWING  COMMANDS  DEFINE  THE  TYPE  OF  LOADING  TO  BE  APPLIED  TO  THE 

*  MODEL.  THESE  COMMANDS  ARE  USED  IN  CONJUNCTION  WITH  THE  TIMEFUNCTION 

*  COMMANDS  DISCUSSED  PREVIOUSLY.  IT  IS  IMPORTANT  TO  NOTE  THAT  IF  THE 

*  APPLIED  LOADING  IS  STATIC  (I.E.  TIME  INDEPENDENT)  THEN  THE  TIME  VALUES 

*  IN  THE  TIME  FUNCTION  COMMAND  REPRESENT  LOAD  STEPS, 

LOADS  TEMPERATURE  TREF=530 

1  1.  0  1  TO  280  1.  0  1 

* 

LOADS  DISPLACEMENT 
3  3  1.0  2 

* 

*  THE  FOLLOWING  COMMANDS  ALLOW  THE  INPUT  NODE  AND  ELEMENT  CONFIGURATION 

*  FOR  THE  MODEL  TO  BE  DISPLAYED  ON  THE  SCREEN  AS  LONG  AS  THE  "WORKSTATION" 

*  COMMAND  IS  COMNE.NTED  OUT. 

FRAME  HEADING=UPPER  XSF=1. 0  YSF=0. 5  XFMAX=22.  5  YFMAX=17.  0 
MESH  N0DES=11  ELEMENT=1 

* 

ADINA 

E.ND 


■# 
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APPENDIX  B:  DATA  REDUCTION  PROGRAMS 


The  following  programs  were  used  to  reduce  the  ADINA  output  file  to  data  points 
to  be  plotted  on  stress-versus-strain  curves. 


C  THIS  PROGRA.M  READS  THE  ADINA  OUTPUT  FILE  AND  COMBINES  THE  STRESS  VALUES 
C  WITH  THE  ELEMENTAL  VOLUMES  TO  PRODUCE  A  VOLUME -AVERAGED  VALUE  OF  AXIAL 
C  STRESS. 

C 

PROGRAM  TEST 
CALL  INIT 
CALL  READ 
STOP 
END 
C 

C  THIS  SUBROUTINE  ZEROS  OUT  ARRAYS  AND  READS  THE  VOLUME  DATA  FILE. 

C  THE  PARA.METER  VARIABLES  ARE:  NUMBER  OF  ELEME.NTS  (NEL),  NUMBER  OF  LOAD 
C  STEPS  (LC)  AND  NUMBER  OF  ELEMENTS  IN  THE  FIBER  (NFBR). 

C 

SUBROUTINE  INIT 

PARAMETER  ( NEL=25  2 , LC=6 , NFBR=1 14 ) 

REAL  SZ( 1:  NEL, 1:  4) ,SZBAR( 1:  LC) ,V( 1:  NEL) ,FBR( 1:  NFBR) 

COMMON  //  SZ,SZBAR,V,FBR 
C 

DO  10  1=1, NEL 
V(I)=0.  0 
DO  15  J=l,4 
SZ(I,J)=0. 0 
15  CONTINUE 

10  CONTINUE 

C 

DO  20  K=1,LC 
SZBAR(K)=0.  0 
FBR(K)=0.  0 
20  CONTINUE 

C 

OPEN  (01,STATUS='0LD’ ) 

DO  25  L=1,NEL 

READ  (01,*,END=50)  V(L) 

25  CONTINUE 

C 

50  PRINT  *,  'EOF  VOL.  DAT' 

C 

RETURN 

END 

C 

C  THIS  SUBROUTINE  READS  THE  FORMATTED  ADINA  OUTPUT  FILE,  AFTER  THE  HEADERS 
C  HAVE  BEEN  REMOVED,  ANT  CALCULATES  THE  AVERAGE  Z  STRESS  IN  EACH  ELEMENT 
C  BY  AVERAGING  THE  FOUR  STRESS  VALUES  AT  THE  GAUSS  POINTS.  THIS  AVERAGE 


C  STRESS  IS  MULTIPLIED  BY  ITS  ELEMENTAL  VOLUME  AND  SUMMED  OVER  THE  ENTIRE 
C  VOLUME  OF  THE  MODEL.  THE  OLTPUT  IS  TWO  FILES:  ONE  CONTAINING  THE 
C  COMPOSITE  STRESS  VALUES  AT  THE  GIVEN  LOAD  LEVELS  AND  ONE  CONTAINING 
C  THE  STRESS  LEVELS  IN  THE  FIBER  AT  THE  SAMxE  LOAD  LEVELS. 

C 

SUBROLTINT  READ 

PARAMETER  (NEL=252 ,LC=6 ,V0L=4.  97,NFBR=114) 

REAL  SZ(  1:  NEL,  1;  4)  ,SZBAR(  1;  LC)  ,V(  1:  NEL)  ,FBR(  1:  NTBR) 

CHARACTER*8  STATE 
COMMON  //  SZ,SZBAR,V,FBR 
C 

OPEN  (02, STATUS= ' OLD ' ) 

OPEN  (03,FILE='V30.  DAT’ , STATUS= ’ NTW ' ) 

OPEN  (04,FILE='V30  FBR.  DAT' ,STATUS='NEV‘ ) 

C 

DO  20  K=1,LC 
SUM=0. 0 
DO  25  L=1,NEL 
READ  (02,*)  lEL 

READ  (02,21)  IN'R, INS, STATE, SX,SY,SZ(L,1),SYZ,SVM,SMAX 
READ  (02,*)  SYS.SMIN 
READ  (02,*) 

READ  (02,21)  INR,INS,STATE,SX,SY,SZ(L,2),SYZ,SVM,SMAX 
READ  (02,*)  SYS,SMIN 
READ  (02,*) 

READ  (02,21)  INR,INS,STATE,SX,SY,SZ(L,3),SYZ,SVM,SMAX 
READ  (02,*)  SYS,SMIN 
READ  (02,*) 

READ  (02,21)  INR,INS,STATE,SX,SY,SZ(L,4),SYZ,SVM,SMAX 
READ  (02,*)  SYS.SMIN 
READ  (02,*,END=50) 

READ  (02,*,END=50) 

21  FORMAT  (14X,I1,3X,I1,5X,A8,5X,4(E12.5,2X),5X,2(E12.5,4X)) 

C 

SBAR=((SZ(L,l)+SZ(L,2)+SZ(L,3)+SZ(L,4))/4.  0)*V(L) 

C 

SUM=SUM+SBAR 

C 

IF  (L.  EQ.NTBR)  THEN 
FBR(K)=SUM/VOL 
END  IF 

25  CONTINUE 

SZBAR(K)=SUM/VOL 
20  CONTINUE 

C 

50  PRINT  *,  'EOF  DETECTED' 

C 

DO  35  N=1,LC 

WRITE  (03,*)  SZBAR(N) 

WRITE  (04,*)  FBR(N) 

35  CONTINUE 

RETURN 
END 


72 


> 


I 


¥ 


C  THIS  PROGRAM  CALCULATES  ELEMENTAL  VOLUMES  BASED  ON  THE  GEOMETRY  OF  THE 
C  FINITE -ELEMENT  MODEL.  THE  REQUIRED  INPUTS  ARE: 

C  Ni::  NUMBER  OF  ELEMENTS  IN  THE  MODEL  IN  THE  HORIZ.  DIRECTION. 

C  NT:  NUMBER  OF  ELEMENTS  IN  THE  MODEL  IN  THE  VERT.  DIRECTION. 

C  NFX:  NUMBER  OF  FIBER  ELEMENTS  IN  THE  HORIZ.  DIRECTION. 

C  NFY:  NUMBER  OF  FIBER  ELEMENTS  IN  THE  VERT.  DIRECTION. 

C  NTM:  NUMBER  OF  ELEMENTS  IN  THE  VERT.  DIRECTION  FROM  THE 

C  MODEL  MIDPOINT  TO  THE  FIBER  END. 

C  NEL:  NUMBER  OF  TOTAL  ELEMENTS  IN  THE  MODEL  (NX  x  NT). 

C  XF:  RADIUS  OF  FIBER  DIVIDED  BY  NFX. 

C  XM:  (RADIUS  OF  CELL  -  RADIUS  OF  FIBER)  DIVIDED  BY  (NX-NFX). 

C  YFl:  (LENGTH  OF  FIBER  -  0.  5*LENGTH  OF  CELL)  DIVIDED  BY  NFM. 

C  YF2:  0.  5*LENGTH  OF  CELL  DIVIDED  BY  (NTY-NFM). 

C  YM:  (LENGTH  OF  CELL  -  LENGTH  OF  FIBER)  DIVIDED  BY  (NY-NFY) 

C 

PROGRAM  CVOL 

PARAMETER  (NX=9 ,NT=27 ,NTX=6 ,NFY=18 ,NFM=4 ,NEL=243 , 

&  XF=.  1666, X.M=.  1633, YF1=.  19,YF2=.  16,YN=.  1644) 

REAL  V0L(1:NEL),V(1:NX,1:NY) 

C 

C  INITIALIZE  THE  MATRICES 
DO  10  1=1, NX 
DO  10  J=1,NY 
10  V(I,J)=0.0 

C 

DO  15  1=1, NEL 
15  VOL(I)=0. 0 

C 

C  CONSTRUCT  THE  V  MATRIX  AND  MAP  IT  INTO  THE  VOL  MATRIX. 

C  TOP  FIBER  GROUP 

NTD=1 


DO  20  L1=NTY,NFY-NFM+1,-1 
DO  21  K1=NFX,1,-1 

V(K1,L1)=  ((XF*K1)**2  -  (XF*(Kl-l))**2)*YFl/2. 
V0L(NTD)=V(K1,L1) 

NTD=NTD+1 
21  CONTINUE 

20  CONTINUE 

C  ***  BOTTOM  FIBER  GROUP  *** 

DO  200  L10=NFY-NFM,1,-1 
DO  210  K10=NFX,1,-1 

V(K10,L10)=  ((XF*K10)**2  -  (XF*(K10-1) )**2)*YF2/2. 
VOL(NID)=V(K10,L10) 

NID=NTD+1 
210  CONTINUE 
200  CONTINUE 

C  ***  TOP  RIGHT  MATRIX  ELEMENT  GROUP  *** 

DO  22  L2=NT,NTY+1,-1 
DO  23  K2=NX,NFX+1,-1 
V(K2,L2)=((XM*(K2-NFX)+XF*NTX)**2- 
&  ( XM*( K2 - 1 -NFX) +XF*NFX ) **2 ) *YM/ 2 . 

V0L(NID)=V(K2,L2) 

NID=NID+1 


23  CONTINUE 
22  CONTINUE 

C  ***  TOP  LEFT  MATRIX  ELE.ME.NT  GROUP  *** 
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DO  24  L3=NT,NTY+1,-1 
DO  25  K3=NFX,1,-1 

V(K3,L3)=  ((XF*K3)**2  -  (XF(K3-1) )**2)*YM/2. 

V0L(NID)=V(K3,L3) 

NID=NID+1 

25  CONTINUE 

24  CONTINUE 

C  ***  MIDDLE  MATRIX  ELEMENT  GROUP  *** 

DO  26  L4=NFY,NFY-NFM+1,-1 
DO  27  K4=NX,NFX+1,-1 
V(K4,L4)=  ((XM*(K4-NTX)+XF*NFX)**2  - 
&  ( XM* ( K4 - 1 - NFX ) +XF*NFX ) **2 ) *YF 1 / 2 . 

V0L(NID)=V(K4,L4) 

NID=NID+1  J 

27  CONTINUE 

26  CONTINUE 

C  ***  BOTTOM  MATRIX  ELEMENT  GROUP  *** 

DO  260  L40=NFY-NFM,1,-1 
DO  270  K40=NX,NFX+1,-1 
V(K40,L40)=  ((XM*(K40-NTX)+XF*NFX)**2  - 
&  ( XM* ( K40 - 1 -NFX ) +XF*NFX)**2 ) YF2/ 2 . 

VOL(NID)=V(K40.L40) 

NID=NTD+1 
270  CONTINUE 

260  CONTINUE 

C  ***  CHECK  TO  ENSURE  THE  SUM  OF  THE  ELEMENTAL  VOLUMES  MATCH  THE  *** 

C  ***  CALCULATED  MODEL  VOLUME  AND  WRITE  THE  OUTPUT  TO  A  FILE.  *** 

VNUM=  ( ( ( XF*NFX ) +( XM* ( NX -  NTX ) ) ) **2 ) * ( ( YF 1*NFM )  + 

&  ( YF  2* ( NFY - NFM ) ) + ( YM* ( NY -NFY ) ) ) / 2 . 

C 

VCALC=d. 0 

DO  30  N=1,NEL 

VCALC=VCALC+VOL(N) 

30  CONTINUE 

PRINT  *,  'VNUM="  VNUTl 
PRI.NT  *,  'VCALC=' ,VCALC 
DO  35  M=1,NEL 
35  WRITE  (10,*)  VOL(M) 

C 

STOP 

END 


C  THIS  IS  A  MODIFICATION  TO  SUBROUTINE  READ,  WHICH 
C  CALCULATES  STRESS  LEVELS  ALONG  THE  FIBER  HALF-LENGTH, 
C  AT  GIVEN  LOAD  LEVELS. 

C 

SUBROUTINE  READ 

PARAMETER  ( NEL-52 , LC-1 , TVOL-3 . 0 , YFBR-1 3 . 0 , XFBR-4 . 0 ) 
REAL  SZ( 1 :NEL, 1 :4) ,SZBAR( 1 :LC) ,V( 1 :NEL) 

REAL  TTFBR(1:LC) ,STFBR(1;YFBR) 

CHARACTER*8  STATE 

COMMON  //  SZ,SZBAR,V,TTFBR,STFBR 
C 

OPEN  ( 02, STATUS-' OLD ' ) 

OPEN  (03, FILE-' FIBER.DAT' , STATUS- ' NEW' ) 

OPEN  (04,FILE-'FIBER  STRESS . DAT ', STATUS- ' NEW ' ) 


74 


I 


tf 


c 

DO  15  K=1,LC 

OPEN  (01,STATUS-'OLD' ) 

FSUM=0.0 

C 

DO  20  M=l,yFBR 
SUM-0 . 0 
VOL-0 . 0 
C 

DO  25  L-1,XFBR 
READ  (02,*)  lEL 

READ  (02,21)  INR, INS, STATE, SX, SY, SZ(L, 1) ,SyZ,SVM,SMAX 
READ  (02,*)  SYS,SMIN 
READ  (02,*) 

READ  (02,21)  INR,INS,STATE,SX,Sy,SZ(L,2) ,SyZ,SVM,SMAX 
READ  (02,*)  SYS,SMIN 
READ  (02,*) 

READ  (02,21)  INR,INS,STATE,SX,Sy,SZ(L,3),SYZ,SVM,SMAX 
READ  (02,*)  SYS,SMIN 
READ  (02,*) 

READ  (02,21)  INR,INS,STATE,SX,SY,SZ(L,4) ,SyZ,SVM,SMAX 
READ  (02,*)  SYS,SMIN 
READ  (02,*) 

READ  (02,*) 

21  FORMAT  ( 14X, II, 3X, II , 5X,A8,5X,4(E12. 5, 2X) , 5X,2(E12 .5, 4X) ) 
READ  ( 01 ,*, END-50 )  V(L) 

SBAR-( (SZ(L,1 )+SZ(L,2)+SZ(L,3)+SZ(L,4) )/4.0)*V(L) 

SUM-SUM+SBAR 

VOL-VOL  +  V( L  ) 

25  CONTINUE 
C 

STFBR(M)-SUM/VOL 

FSUM-FSUM+SUM 

WRITE  (03,*)  STFBR(M) 

20  CONTINUE 
C 

CLOSE  (01) 

TTFBR( K )=FSUM/TVOL 

WRITE  (03,*)  'LOAD  STEP-',K 

WRITE  (04,*)  TTFBR(K) 

15  CONTINUE 
C 

CLOSE  (02) 

50  PRINT*,  'EOF  DETECTED' 

RETURN 

END 
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